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Abstract 

Quasi-normal-eigenvalue optimization is studied under constraints 6 i(x) < B{x) < b 2 {x) on 
structure functions B of 2-side open optical or mechanical resonators. We prove existence of 
various optimizers and provide an example when different structures generate the same optimal 
quasi-(normal-)eigenvalue. To show that quasi-eigenvalues locally optimal in various senses are 
in the spectrum of the bang-bang eigenproblem y" = —uj‘^y[bi -|- (62 — where 

indicator function of C+, we obtain a variational characterization of in terms 
of quasi-eigenvalue perturbations. To address the minimization of the decay rate |Ima;|, we 
study the bang-bang equation and explain how it excludes an unknown optimal B from the 
optimization process. Computing one of minimal decay structures for 1-side open settings, 
we show that it resembles gradually size-modulated 1-D stack cavities introduced recently in 
Optical Engineering. In 2-side open symmetric settings, our example has an additional centered 
defect. Nonexistence of global decay rate minimizers is discussed. 


MSC-classes: 49R05, 47J10, 35B34, 49M05, 46N10, 90C29, 47A55, 78M50, 70J10, 34L15 

Keywords: high-Q, photonic crystal, quality-factor, Pareto optimal structural design, resonance 
perturbation 


1 



2 


Contents 


1 Introduction 


2 

3 


4 


5 


6 


7 

8 

9 

10 


11 

12 


A 


Basic settings, quasi-eigenvalues and admissible sets 

Properties of resonances and related maps 

3.1 Characteristic determinant F and definition of multiplicity . . . 

3.2 Non-simple resonances, the case of nonnegative B, and examples 

Perturbations and weak continuity of resonances 

4.1 Local weak continuity. 

4.2 Resonances’ perturbations and directional derivatives of F ... 

Variational characterization of the nonlinear spectrum 

5.1 The turning interval and rotational properties of modes. 

5.2 Proof of Theorem 15. II . 


Local extremizers and optimization for a fixed frequency 

6.1 Definitions of various optimizers. 

6.2 Admissible frequencies and examples of admissible families 

6.3 Existence of optimizers and properties of S[A] . . 

Non-uniqueness and local minima for small a 
Local minima in the case of small contrast 

Locally extremal uj are nonlinear eigenvalues 

9.1 Local boundary of S[A] and the proof of Theorem 

5-D reduction with the use of bang-bang equations 

10.1 Finding of optimizers from the nonlinear spectrum. 

10.2 Nonlinear eigenvalues that are not of minimal decay. 

10.3 Nonlinear eigenvalues as zeros of a function of three variables 

A numerical experiment for constant side constraints 
Discussion 

12.1 Alternating almost periodic structures with a centered defect. 

12.2 Global minimization of decay rate. 


9.1 


Appendix 

A.l Analyticity, the proof of Lemma 3^ 

A.2 Directional derivatives, the proof of Lemma [4^ 

A.3 Admissible high frequencies, the proof of Proposition 6.5 
A.4 Solvability of bang-bang equations, the proof of Theorem |10.4 


111 


m 

m 

m 

[9] 

IQ] 

m 

m 

12 ] 

m 

IS] 

H] 

13 

n 

m 

[ 21 ] 

[23] 

[25] 

[25] 

[26] 

m 

[23 

[23 


m 

Eo] 

m 


[32] 


m 

m 


El 




























3 


1 Introduction 


The paper is devoted to the analytical and nnmerical stndy of two connected qnestions: nonlinear 
eigenproblems with a bang-bang term and optimization of qnasi-(normal-)eigenvalnes in optical and 
mechanical resonators. 

A leaky cavity (or resonator) is a region of space within which the electromagnetic (or aconstic) 
held is well conhned, bnt not completely conhned. Cavities snpporting eigenmodes with high qnality 
factor (high-Q) are needed for a nnmber of applications inclnding Cavity Qnantnm Electrodynamics, 
Optical Engineering, and Microscopy (see e.g. [13 |3H ESI SO])- The recent progress in fabrication 
of small size optical resonators [H ESI EH] attracted considerable interest to nnmerical [IHIESIEIEII 
EHIEI] and analytical [SHEOlES] aspects of resonance optimization. 

Optical qnasi-eigenvalnes (or resonances) can be mathematically dehned as eigenvalnes u cor¬ 
responding to ontgoing solntions of the time-harmonic Maxwell eqnation [CT EH]. Becanse of the 
leakage, resonant eigenoscillations are expected to decay exponentially in time. The modulns of the 
imaginary part | Ima;| of a qnasi-eigenvalne corresponds to the decay rate of the standing wave, the 
real part Rea; to the frequency of oscillations. In Engineering stndies, the conhnement of energy for 
a particnlar resonant eigenmode is measnred by the quality factor Q{u) = 

Since light is essentially difficult to localize it is hard to realize small-sized optical cavities with 
strong light conhnement. One of the reasons of active involvement of Numerics in this held is that the 
initial progress with fabrication of high-Q cavities based on 2-D photonic crystals was achieved with 
the help of computer simulations of eigenmodes for perspective versions of design [1]. The interest 
to multilayer structures with 1-D geometry has soon returned [331 EH EHl El EU, partially because 
of cheaper fabrication process, partially because computations in this case are simpler. For ideal¬ 
ized 1-D cavities with inhnite layers and under the assumption of normally passing electromagnetic 
(EM) waves, the Maxwell system can be reduced to a 1-1-1 wave equation of nonhomogeneous string 
B{x)dfu{x,t) = dlu{x,t) (see e.g. |l3l SO]). The analogy with a string is widely used in Optical 
Engineering [271 EH] (and possibly goes back to Lord Rayleigh [TTl Chapter 4]). 

Quasi-eigenvalues associated with B{x)d^u{x,t) = dlu{x,t) can be dehned via the eigenproblem 
consisting of the time-harmonic equation 

y”{x) = —upB{x)y{x) a.e. on oi < a; < 02 , (1.1) 


equipped with radiation boundary conditions. In this 1-D case, the radiation conditions take the 
a;-dependent form of local conditions at the endpoints 


y\ai) 

—iuj 

y'{a2) 

iUJl'2 


I'lViai), 

y{cL2) (see Section ^for details). 


( 1 . 2 ) 

(1.3) 


The following properties of the set E(i?) of quasi-eigenvalues associated with the coefficient B{-) 
are important for the present paper: E(i?) is a subset of the lower complex half-plane C_ = {z E 
C : Imz < 0} and is symmetric with respect to (w.r.t.) the imaginary axis iM, quasi-eigenvalues are 
isolated and of hnite algebraic multiplicity, 00 is their only possible accumulation point (see Section 
l^and, for one-side open cases, [HU Ell El |l2l [36]L 

From mathematical point of view, estimates on resonances associated with various wave equations 
have being studied in Mathematical Physics at least since 1970s [2H| and is still an active area of 
research (see e.g. [2] and references therein). Optimization of resonances may be seen as an attempt 
to obtain sharp estimates of this kind. This point of view and the study of resonances associated 
with random Schrodinger operators were initial sources of the interest to the problem HEllIll. It 
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seems that, up to now, sharp estimates have been obtained only in low frequency regions for a few 
1-D models involving total mass type constraints [201 ES] (on a somewhat different problem involving 
damping, see [ 6 ] and the discussion in 0 )- Among other variational problems for eigenvalues, the 
optimization of resonances can be classified as the nonselfadjoint spectral optimization. It contains 
essentially new effects and difficulties in comparison with selfadjoint cases (see [3 El H EB E3]). 

Since quasi-eigenvalues u are complex numbers, one of the ways to formulate related variational 
problems is to interpret u as an M^-vector and to take the point of view of the optimization theory 
for vector-valued cost functions (Pareto optimization) [23]. The approach of [23] is a development of 
that of [El]) where the engineering problem of the high-Q design for one-side open multilayer cavities 
was considered analytically. To provide a simple and simultaneously rigorous formulation, the papers 
[201 EB restate the problem of high-Q design in terms of the minimization of the decay rate | Ima;| 
for quasi-eigenvalues with a hxed frequency Reca = a. Then the set of all quasi-eigenvalues with 
minimal possible decay for their particular frequencies forms a Pareto optimal frontier (see Fig. [Tlin 
Section § and [231 Fig. 1]). 

The quasi-eigenvalues are supposed to be generated by cavities with structures belonging to a 
certain admissible family. In the present paper, the resonator’s structure is described by an L^- 
function B{x) on [ 01 , 02 ] (—00 < oi < 02 < + 00 ), which represents spatially varying dielectric 
permittivities in the case of optical cavity, or the linear density in Mechanics settings. Boundary 
conditions at the endpoints 01 ^ 2 , which describe the leakage of energy due to radiation or damping, 
include two parameters z/i, z /2 G [0, -|-oo]. These parameters are assumed to be fixed, and so, they do 
not participate in the optimization process. The admissible family of the structures A (over which 
the optimization is performed) is defined by the functional side constraints 61 (x) < B{x) < b 2 {x) a.e. 
on ( 01 , 02 ). This allows us to consider an important for applications situation when the resonator 
itself is only a part of a more complex device [32l |29] . That is, sometimes only some parts of the 
system are suitable for modihcations to achieve better resonant properties. The part of the device 
that is not subjected to optimization corresponds to the intervals where &i(x) = 62 ( 2 ^)- 

The connection between quasi-eigenvalues of minimal decay and nonlinear eigenproblems was 
established recently in [211 Eg for the case of 1-D models of one-side open resonators and constant side 
constraints hi and 62 - R was proved that the set of Pareto optimal quasi-eigenvalues of minimal decay 
belongs to the spectrum associated with the equation y'\x) = —u‘^y{x) bi -\- [62 — ( 2 /^(^)) 


The nonlinearity contains the function of bang-bang type B{y){x) = bi-\- (62 — bi)x^^{y'^) 


taken 


of the square of complex-valued solution y. Here and below Xc^(') is the indicator (characteristic) 


function of the open upper complex half-plane C+ = {z : Imz > 0}, i.e.. 


^C4 


= 1 if z G C+, and 


= 0 if z G C_ = {z : Imz < 0}. 


So B{y){-) equals 6 i(-) on the x-intervals where \my‘^[x) < 0, and B{y){-) = h 2 {-) for x-intervals where 
\my‘^{x) > 0. This resembles to some extend bang-bang equations arising in the optimal control 
theory when the control is restricted by a lower and an upper bound. However, the interpretation of 
the above equation from the point of view of control is not found yet. 

Since the nonlinear equation derived in [21] excludes infinite-dimensional unknown structure 
R, it may essentially simplify computation of optimal resonances in comparison with step-by-step 
adjustment procedures for B based on computation of electromagnetic held on each step. 

Development and analytic justihcation of a numerical procedure based on the bang-bang eigen- 
problem is one of the goals of our research. However, the implementation of this new computational 
approach requires better understanding of the shapes of the nonlinear spectrum and of the Pareto 


optimal set, as well as their projections to the real line (see Section 6.2). Our method also requires 


the study of properties of solutions of the bang-bang equation (1.4), which are considered in Section 


10.3 and Appendix A.4 
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We show that not only qnasi-eigenvalnes of minimal decay, bnt also qnasi-eigenvalues locally 
extremal in a mnch wider sense belong to the set of eigenvalnes associated with the nonlinear 
eqnation 


y"{x) = -u^y bi{x) + [b 2 {x) - bi{x)]x^^{y^{x)) 


Oi < X < 02 


(1.4) 


(which inclndes the fnnctional side constraints &i, 2 (') and is snpposed to be eqnipped with the ra¬ 
diation boundary conditions). In particular, we show that nonlinear spectrum contains the 
nonimaginary part bdS[A] \ iM of the boundary bdS[A] of the set of admissible quasi-eigenvalues 


E[A| := U E(B) 


see 


Sections and . 


seA 


To achieve this, we provide a variational characterization of the part of lying in C \ iM (Theorem 
5.1). We also show that bdS[A] and so are, at least in some cases, essentially ‘greater’ than the 
set of quasi-eigenvalues of minimal decay (see Sections and 11 and Example 10.3). 

The quasi-eigenvalue of minimal decay for a particular frequency a can be found as the closest 


to M point of the intersection of the nonlinear spectrum with the line a - 1 - iM (see Section 10). Then, 
corresponding optimal structures B can also be found from the bang-bang eigenproblem. In the end 
of the paper, we support our study by a numerical experiment (see Fig. [2| (a) and Tables 1-2 in 


Section 11). 


To prove rigorously that points of the nonlinear spectrum are not necessarily quasi-eigenvalues of 
minimal decay, the case of small dielectric contrast is considered in Section In Optics literature the 
study of small contrast is often serves as a base of intuitive understanding of qualitative effects (see 
e.g. [m Chapter 4]). In the optimization settings, small contrast means that the constraints bi and 
62 are close to each other. To measure the corresponding distance we use the L^-norm ||52 since 

this choice is flexible enough to cover both the small changes in permittivity function (when || 62 ~^i||oo 
is small) and fluctuations in the widths and positions of layers. This gives a connection with models 
involving uncertainties or random deviations in the resonator’s structure B{x). Indeed, for stochastic 
models, all possible positions of random resonances (the supports of probabilistic distributions) can 
be localized by the set of admissible resonances S[A] of a properly posed optimization problem. 
In practical situations uncertainties or fluctuations in the cavity’s structure may model unintended 
deviations in fabrication process [16], incorporation of biomolecular thin Elms in X-ray multi-layer 
cavities [35], and various Physics effects affecting EM properties of materials [T 6 ] . 

Pure imaginary quasi-eigenvalues (i.e., the case u G iM) correspond to the over-damping and 
critical damping effects (in Mechanics settings). For them, optimization partially resembles that for 
selfadjoint problems [SIEIEIEI]. The pure imaginary case is considered separately in Section It 
provides, as a by-product, an example of two different optimal structures associated with one quasi¬ 
eigenvalue of minimal decay. This partially answers the uniqueness of optimizer question discussed 
in [IllEniE3|. 

The application of steepest ascent numerical procedures to the search of quasi-eigenvalues with 
locally minimal decay rate (or locally maximal quality factor) attracted a considerable attention 


dsiiisiiini. It seems that Corollary |7.3| (i) and Theorem | 8 .1 1 provide first rigorous proofs of existence 
of such local minimizers. (The proof of Theorem 8.1 can be easily transformed to consider local 
maximizers for the Q-factor). The existence of global minimizer for decay rate (over all G A and 
oj G C_ such that u G E(i?)) is questionable. Our opinion is that it is difficult to expect such 


existence in most of practical situations, see the discussion in Section 12.2 


The discussion section (Section 12) also compares the results of the numerical experiment of 
Section IT with optimal design suggestions of engineering and numerical optimization papers. Our 
conclusion is that the gradually size-modulated 1-D stack designs of janiEniE] are reasonable, but 
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for symmetric resonators, it is possible that an additional defect in the center may improve the 
confinement of energy. 

The results of the paper were partially reported at the International Congress of Mathematicians 
(Seoul ICM 2014) and in several seminar and workshop talks at the University of Liibeck, the Institute 
of Mathematics (at Kiev), B. Verkin Institute for Low Temperature Physics and Engineering, and 
Pidstryhach Institute for Applied Problems of Mechanics and Mathematics. 

Notation. By [x] and [xj the ceiling and floor functions are denoted, i.e., the smallest integer 
not less than x and the greatest integer not greater than x, respectively. 

The following sets of real and complex numbers are used: open half-lines = {x G M : ±x > 0}, 
and open discs ©^((C) := {z E C : \z — C\ < e} with the center at ( and radius e. The Lebesgue 
measure of a set S' C M is denoted meas S. 

L^(g)( 01 , 02 ) are the Lebesgue spaces of complex- (resp., real-) valued functions and 

lU''’^[oi, 02] := {y E L^(oi, 02) : d^y E L^(oi, 02), I < j < k} 

are Sobolev spaces. The corresponding standard norms are denoted by || ■ ||p and || ■ \\wk,p. The space 
of continuous complex valued functions with the uniform norm is denoted by C[ai, 02 ]. 

Under the support suppi? of 5 G L^(oi, 02 ) we understood the topological support of the ab¬ 
solutely continuous measure |i?(x)|dx on [ 01 , 02 ], i.e., suppi? is the smallest closed set S such that 
B(x) = 0 a.e. on [oi, 02 ] \ S. 

For basic definitions of convex analysis we refer to [39]. Let S' be a subset of a linear space U over 
C (including the case U = C). For Uq E U and z E zS + Uq := {zu + Uq : u E S'}. The convex 
cone generated by S (the set of all nonnegative linear combinations of elements of S') is denoted by 
cone S'. Open balls in a normed space U are denoted by 


Be(Mo) = ®e(Mo; U) -.= {u E U : 


\u 


< e}- 


The closure of a set S (in the norm topology) is denoted by S, the boundary of S by bdS'. 
For a function / dehned on S', /[S'] is the image of S'. By d^f, d^f, etc., we denote (ordi¬ 
nary or partial) derivatives w.r.t. x, z, etc. The function /+ is ‘the positive part’ of f, i.e., 
/(x) , if /(x) > 0 _ 

0 , if /(x) < 0 ’ 

X G S', and Xs'(x) = 0 when x ^ S'. 


f+{x) := 


Xs is an indicator function of a set S', i.e., Xs{x) = 1 when 


2 Basic settings, quasi-eigenvalues and admissible sets 


The boundary conditions (1.2) and ( |1.3[ ) at 01^2 £ 1^ involve the (extended) constants ui^ 2 , which 
throughout the paper are assumed to satisfy 


1^1 G [0,-Fcx)), z /2 G (0,-Fcx)], |z/i||l/i/ 2 | 7 ^ 0, and z/i < z/ 2 . 


( 2 , 1 ) 


When damping coefficients 1^12 are fixed, we can define the set of quasi-eigenvalues F(i?) associated 
with the structure B as the set of a; G C \ {0} such that eigenproblem (1.1), (1.2), (1.3) has a 
nontrivial solution y E lU^’^[ 01 , 02 ] (i.e., a solution that is not identically zero). This solution y is 
called a (quasi-normal) mode. Several other names for u are used, sometimes in slightly different 
settings: scattering poles [28], dissipation frequencies [21125], resonances [H HH] , quasi-normal levels 
(in the Quantum Physics literature). 

The function B E L^{ai, 02 ) in (1.1) describes the structure of a non-homogeneous medium vary¬ 


ing in x-direction. In the case of optical cavity, B[xq) corresponds to permittivity of the transparent 
dielectric material of the layer with the x-coordinate equal to xq (if the light speed in vacuum is 
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normalized to be 1 or if one works with the complex wave number instead of the complex angular 
frequency u). In Mechanics models involving the equation of nonhomogeneous string, the function B 
is the varying linear density of the string (if the tension equals 1). So, in these models, B is supposed 
to be a.e. nonnegative, and we will keep this assumption in the context of optimization problems. 
However, in Section the dehnition of quasi-eigenvalues will be extended to complex-valued coeffi¬ 
cient B for the needs of perturbation theory of Section 

When Pj is in (0, -|-oo), the corresponding boundary condition is dissipative and describes either 
linear damping at Oj or the radiation of waves into outer regions (see the examples in Section 6.2 
for details). Our settings for a 1-D resonator is a slightly generalized version of 


We allow 

B to be equal to 0 on certain sets (the massless string approximation). We also include the cases 
when one of the boundary conditions is conservative (and so the resonator is only one-side open). 


In particular, when P 2 = +C) 0 , we suppose that 1 / 1^2 = 0, and then (1.3) turns into the Dirichlet 


condition 7 /( 02 ) = 0. Note that the assumption Pi < P 2 does not restrict the generality since, in the 
case Pi> P 2 i one can use the change of variable x = —x. 

While most of mathematical studies were devoted to one-side open or symmetric resonators 


(mathematically, these two cases are almost equivalent, see Example 6.4), contemporary Optics ap¬ 
plications often involve the two-side open case [33l 021 [26], jS] [31]. In the present paper, we consider 
one- and two-side open resonators in a unihed way. Besides the importance for engineering appli¬ 
cations, the study of two-side open resonators brings new mathematical effects. In particular, the 
spatial phase axgy{x) of the complex-valued mode y is not monotone in contrast to one-side open 
case. This effect, and the fact that the notion of switch points (see 121 ) looses its natural sense when 
the constraint functions 61^2 are allowed to be equal on a set of positive measure, lead in Section]^ to 
more detailed study of rotational properties of resonant modes and to more essential use of Convex 
Analysis in comparison with the technique of [ 21 ] . 

We consider optimization over the following family of admissible structures 


A := {B{x) G L^{ai,a 2 ) : bi{x) < B{x) < 62 ( 2 ^) a.e. on ( 01 , 02 )}, 
where 61^2 are certain Lebesgue integrable functions dehned on ( 01 , 02 ) such that 

0 < bi(x) < 62 ( 37 ) on ( 01 , 02 ) and 
measA > 0, where A = {x G ( 01 , 02 ) : bi(x) < b 2 (x)}. 


( 2 , 2 ) 

(2.3) 

(2.4) 


A complex number uj is called an admissible quasi-eigenvalue if it belongs to the set S[A] := 
Tj{B). Keeping in mind applications, we will pay main attention to the problem of minimization 
of decay rate of an individual resonance u generated by a certain B & A. However, from mathematical 
point of view it is more convenient to consider u itself as a cost function with the values in and 
to study extremal in various senses quasi-eigenvalues . One of such sets is the boundary bdE[A] of 
the set of admissible quasi-eigenvalues E[A]. 

Actually, u can not be considered directly as a functional of B even locally due to multiplicity 
and splitting issues (see Sections 3.2| Proposition [4.4[ and the discussion in [23|). However, a rigorous 
approach to u:{B) can be given via the set-valued map B 1 —)■ E(H). Local extrema of such set-valued 
maps were introduced in [23| and will be the key tool in Section In such generalized settings, 
bdE[A] plays a role of a set of generalized Pareto extremizers for the map B 1 —)■ E(i?) over A. 


3 Properties of resonances and related maps 

3.1 Characteristic determinant F and definition of multiplicity 

In this section, an equivalent (but more convenient from the point of view of perturbation theory) 
dehnition of quasi-eigenvalues is given under more general assumption that B G Lc(ai,a 2 ). 










Denote by (p{x) = ip{x,z;B), 'ijj{x) = z] B), and 6{x) = 9{x,z;B) the solutions to dly{x) = 
—z‘^B{x)y{x) satisfying 

(p{ai,z] B) = 9x^(ai, B) = 1, da,(p{ax, z] B) = i>{ai, z] B) = 0, 

6*(ai, z;-B) = 1, dx0{ai^ z] B) = —izvi. (3.1) 

Obviously, 6{x) = y){x) — izi'i'il){x) and 6{x) is a unique solution to the integral equation 

y{x) = 1 — izui(x — ai) — z‘^ / {x — s) y{s) B{s)ds, Oi < x < 02 - (3.2) 

J ai 

Recall that the set of quasi-eigenvalues corresponding to a structure B (in short, quasi-eigenvalues 
of R) is denoted by S(i?). The following lemma is the integral reformulation of the quasi-eigenvalue 
problem. 

Lemma 3.1. A number u E C belongs to S(R) if and only if there exists nontrivial y{x) G O[0,£] 
satisfying equation (3.2) with z = u and the equality 

'■02 


y{a2) + ^ 

^2 T^2 


y{s) B{s)ds = 0. 


(3.3) 


Oi 


Note that in the integral settings there is no need to exclude separately the case ca = 0. If a; = 0, 
then it is easy to see that the problem consisting of (3.2) with z = u and (3.3) has no nontrivial 
solutions. 

Clearly, S(i?) is the set of zeroes of the function F{--, B), where 

) 

6{s, z] B) B{s)ds. 

^2 1^2 Jai 

It is easy to see that for 2 ; 7 ^ 0, 


F{z) = F{z; B) := 0 (a 2 , R) + - - - 

^2 ^2 


(3.4) 


F{z;B) = 9{a2,z;B) + 


idx9{a2,z;B) 


Z 02 


(3.5) 


When 1^2 = 00 , these formulas turns into F{z; B) = 6 ^( 02 , z( B). 

We say that a map G : Ui ^ U 2 between normed spaces Ui ^2 is bounded-to-bounded if the set 
G'[S'] in U 2 is bounded for any bounded S in Ui. 

Lemma 3.2. (i) The functional F{z]B) is analytic on the Banach space C x L^(ai,a 2 ). 

(ii) Let I < p < 00 . The mappings {z,B) 1 —>• ip{-,z;B), {z,B) 1 —>■ 'il!{-,z;B), and {z,B) 1 —)■ 
9{-,z;B) are bounded-to-bounded and analytic from C x £^( 01 , 02 ) to W^’^[ai,a 2 ]. 


The proof is given in Appendix A.l (concerning analytic mappings on Banach spaces, see [37]). 

It is obvious that all modes y corresponding to a; G S(R) are equal to 9{-,u; B) up to a multipli¬ 
cation by a constant. So the geometric multiplicity of any quasi-eigenvalue equals 1. In the following, 
the multiplicity of a quasi-eigenvalue means its algebraic multiplicity. 

Definition 3.1. The multiplicity of a quasi-eigenvalue of B is its multiplicity as a zero of the entire 
function F{-]B). A quasi-eigenvalue is called simple if its multiplicity is 1. The set of non-simple 
quasi-eigenvalues is denoted by Smuit(7?) (non-simple quasi-eigenvalues are often called multiple, or 
degenerate). 

This is essentially the classical M.V. Keldysh definition of multiplicity for eigenvalue problems 
with an eigen-parameter in boundary conditions (note that F coincides with the characteristic de¬ 
terminant of ( |1.1 )-(1.3) up to a constant), see [2ll |25l [T] [36l ISU |23] and references therein. 

Since F(0; R) = 1 -|- ^ > 0, each quasi-eigenvalue has a hnite multiplicity and the set D(R) 
consists of isolated points, which can accumulate only to 00 . Note that D(R) may be empty, see 


Proposition 3.4 
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3.2 Non-simple resonances, the case of nonnegative B, and examples 

There exist triples {ui,U 2 ,B), consisting of numbers Ui 2 satisfying the assumption ( |2.1 ) and a non¬ 
negative function B, that generate non-simple quasi-eigenvalues. This follows from [12] and [SUl 
Theorem 4.1] (with V 2 = +oo). For a slightly different class of quasi-eigenvalue problems, existence 
of degenerate quasi-eigenvalues was proved in [251 Theorem 3.1] (English translation of this theorem 
can be found in izni, see also 1211). and examples were given in [aniEni. 


Lemma 3.3. Let B{x) >0 a.e. and ^ 1^2 satisfy (2.1). 


(i) S(5) C C_, 


Then: 


(a) E(i?) is symmetric w.r.t. the imaginary axis iM, moreover, the multiplicities of symmetric quasi¬ 
eigenvalues coincide. 


Proof, (i) Let u G E(i?) and y be an associated eigenfunction. The energy of the eigen-oscillation 
u{x,t) = e~^‘^^y{x) at the time t is 


r-a2 


'■02 


\dxu\‘^ B\dtu\‘^dx = e where En^, := / \y'\‘^B \u:\^ \y\^dx. 


'Ol 


'01 


Since y is nontrivial, we see that En^ > 0. (Otherwise, y’ = d and, in turn, ?/ = 0). 
Recalling that a; 7 ^ 0 and taking the real part in the identity 


r-a2 


r-a2 


—iu En,, = — io; 


\y \ dx -\-iu {—u )B y y dx = 


'ai 


'Ol 


r*a2 


= —2iRea; 


\y'fdx - IwF iyi\y{ai)f + 


'ai 


\y'i^2)? 

|a;Pz/2 


one gets 

Ini a; 

Note that inequality h'i\y{ai)\‘^ -\- 
(ii) follows from identities 


KU 


En„ 


o|2/(ai)P + 


\y'ia2)? 

|n;Pi/2 


< 0 . 


> 0 follows from ( 2 . 1 ) and y ^0. 


iC; B) = e{x, iC; B) and F(iC; B) = F(iC). 


(3.6) 


□ 


When the medium of a resonator is homogeneous, quasi-eigenvalues can be calculated explicitly. 


Proposition 3.4. Let B{x) = b, where h>d is a constant. 

(i) If = z/i or 6 ^/^ = 1 / 2 , then S( 6 ) = 0. 

(a) lfb>0, 6 ^/^ 7 ^ z/i, and 6 ^/^ 7 ^ z/ 2 , then S(&) = {oi}n{b)}f^_^ with 


LVnib) = 


2b^P{a2 — ai 


In 


l + i^i 


1-iFi 


+ 


_ L ^ [ 0 , 1 ^ 2 ] 

b^B{a 2 -ai) \(n-M/ 2 ), if b^^"^ e {^ 1 , 1 ^ 2 ) 


(3.7) 


where Ki := ^ g 77, G Z. 

t . 1/2 ^2 

(ill) Ifb = d^vi, then S( 6 ) = }■ 

In the cases (ii)-(iii) all quasi-eigenvalues are simple. 
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Proof. Consider the case 6 > 0. The equation F{z) = 0 can be transformed into tan([a 2 — ai\F^‘^z) + 
iKi = 0 and we know that the roots u are in C_. Taking dz one can see that the roots are simple. 
We have Ki = 1 exactly when 6 ^/^ G {z/i, ^ 2 ]- In this case S( 6 ) = 0 . In the case Ki 7 ^ 1, the real 
part of Arctan(—ii^i) is either vrn or 7 r(n + l/2) depending on the sign of \ — Ki. This leads to (3.7). 
Note that 


1 > iCi -v^ 6^'^^ ^ 1 < Ki 6^/^ G (z/i,i/ 2 ), and Ki = 1 G {i^i,z^ 2 }- (3.8) 

The case 6 = 0 can be treated by straightforward computations. □ 

In the particular cases of one-side open resonators (i.e., when z/i = 0 or z /2 = -l-oo) this example 
is well known ISIE]. 


4 Perturbations and weak continuity of resonances 


Let [/ be a Banach space. For a functional G{z] u) that maps C x f/ to C, we denote by 


dG{z] uq) 


(mi) := lim 
du C^o 


G{z]Uq + C^i) - G{z]Uo) 

c 


the directional derivative of G{z]-) along the vector ui G f/ at the point Uq G U. By 


dG(z;uo) 

du 


the 


corresponding Frechet derivative in u is denoted. If the Frechet derivative exists, it is a linear 

functional on U and [>5'] defines the image of a set S G U under this functional. 

Remark 4.1 (Global discontinuity oi B ^ S(i?)). Before the study of continuity properties of quasi¬ 
eigenvalues, let us note that new quasi-eigenvalues can appear from 00 with small variations of B. 
To see this, it is enough to take the homogeneous structure Bq equal to uf or z/| and to consider 


S(S) using (3.7) and taking B = Bq + e with small e G 


This case is very special because, when Bq = izj, the resonator is not separated from outer space 
by any reflecting barrier (see detailed explanations in [7]). We do not know if the above effect of 
global in u discontinuity of the set-valued map B 1 —)■ F(i?) can happen in a vicinity of Bq ^ 1 /^ 2 - 


4.1 Local weak continuity 

Let us fix a countable family of continuous functions that is dense in C'[ai,a 2 ]. This family 

generates a metric pM(dMi,dM 2 ) := ^ 2^(1 + [^^dM )|) ^ complex 

Borel measures on [oi, 02 ]. The weak* topology on any closed ball in M coincides with the topology 
generated by the metric pm- We will use the metric 

p(5i, B 2 ) := pmiBidx, B 2 dx) on Lc(ai, 02 ), 


where Bj dx are absolutely continuous measures corresponding to the functions Bj G L\.. 

In this subsection, is the closed ball Bj:j(0) in Lj.(ai,a 2 ) with radius R> d and centered in 
the origin. 


Lemma 4.1. Let z G C. 

(i) The functional F{z] ■) : B^ C is continuous w.r.t. the metric p. 

(a) The map B 1 —>■ 6 *(-, 2 ;; B) is continuous from the metric space (Bi?, p) to the normed space G[ai, 02 ]. 
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Proof, (i) follows from (ii). Let us prove (ii). Suppose B^,Bj G Mr and p{Bj,B^) —)■ 0 as j —)■ cxd. 
In particular, w*-limi?jda; = B^dx. By Lemma 3.2 (ii), 

the sequence 6j{-) := 6{-,z]Bj) is bounded in W^’^, and so is relatively compact in C[ai,a 2 ]. (4.1) 


Hence there exists a subsequence 9j^ (strongly) convergent in C'[ai,a 2 ] to a certain 9^. Using 
w*-limHjda; = i?*da; and passing to the limit in (3.2) for every x G [ 01 , 02 ], one can show that 


9,i-) = 9{;z-,B,). 

Assume that 9j does not converge to 9.^ 
convergent to 0** 7 ^ 0*. Passing to the limit in (3.2) again, we see that 9., 


contradiction. Thus, 9{-, z; Bj) —)■ 9{-, z; B^) as j —)■ 00 


Then (4.1) imply that there exists a subsequence 9n^. 

{.) = 9{;z-B,) = 9,{-), a 

□ 


The total multiplicity of quasi-eigenvalues of B in a set V G C is the sum of multiplicities of all 
OJ G S(H) n V. 

Proposition 4.2. Let Bq G Lj.(oi,02) and R > ||i?o||i- Let V be an open bounded subset of C such 
that its boundary hdV does not contain quasi-eigenvalues of Bq. Then there exists a neighborhood 
W C Bii; of Bq in the topology of the metric space (Br, p) such that, for any B G W, 


there are no quasi-eigenvalues of B on hdV, 
the total multiplicity of quasi-eigenvalues of B inV coincides with that of Bq. 


Proof. By Lemmas 3.2 and 4.1, the functional F{z-, B) is analytic in z and B and is continuous in the 
second variable on the metric space {Mr,p). This allows us to apply [H Theorem 9.17.4] (which is a 
corollary from Rouche’s theorem) to F{z]B). The equivalence z G S(R) F{z-,B) = 0 completes 
the proof. □ 


4.2 Resonances’ perturbations and directional derivatives of F 

Lemma 4.3. At quasi-eigenvalues u G Ll(i?), the derivative of F w.r.t. z is given by 


dF{u;B) 

dz 


2u 


ra,2 


da;9{a2,uj,B) 


9‘^{s,u; B)B{s)ds + 




’ai 


dJ{a2,uj,B) 


+ 


9{a2,uj, B) 


u 


(4.2) 


and the directional derivatives [ObF^uj, B)] (U) w.r.t. B in the direction V G L“(ai,a 2 ) by 


dF{u, B) 

Fb 


(^) = 


bJ 


da:9{a2,UJ,B) 


6 *^(s, u; B) V (s) ds. 


'ai 


(4.3) 


The calculation of derivatives resembles that of 


and is given in Appendix A .2 


Proposition 4.4. Let u be a quasi-eigenvalue of B G of multiplicity m G N. Then: 

(i) For every direction U G ? there exist e > 0 , 5 > 0 , and functions flj, j = 1,... ,m, continuous 
on [0,(5) such that for ( G [0,(5), all the quasi-eigenvalues of B -\- (V lying in D£(a;) are given by 
{r 2 j(C)}]" taking multiplicities into account. 


(ii) For each of these functions the following asymptotic formula is valid 


n,{0=u, + lK{u,B-,V)(Y/"'+o{CG as(^0, with K{u, B-,V) := - g) ■ (4.4) 

(Hi) In the case K(u, B;V) ^ 0, each branch corresponds to exactly one of functions Ilj. So 

for small enough C > 0 , all m values of functions Uj(C) distinct and simple quasi-eigenvalues of 
B + CV. 
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Proof. For a one-side open resonator, an analogue of the particular case considered in (iii) was 
obtained in m Proposition 3.5]. Statement (i) is well known for resonances of the Schrodinger 
equation and can be obtained using the Weierstrass preparation theorem and the Puiseux series 
theory. We start from these arguments to prove that (4.4) holds always (not only when K{u, B, V) ^ 

!)• 


0 as in 

Consider the entire function Q(z, Cf) := F{z] B -|-CC) of two complex variables z and C,. Then uj is 
an m-fold zero of the function Q{-, 0). By the Weierstrass preparation theorem, in a certain polydisc 
V^{uj) X Vs^{0), 


C) — [{z ~ + hi{C){z — a;)™ ^ -f • ■ ■ -I- hjn{C)]R{z, C), 

where the coefficients hj (the function R{z,()) are analytic in Vsj^{u)) (resp., in V^{uj) x T> 5 ^( 0 )), 

hj(0) = 0, and R{z,C,)l^^ in ^ RsAR)- 

Differentiating Q by C, one gets = c)^Q(a;,0) = h(„(0)i?(a;, 0). On the other side, 

a™F(a;, B) = df^Q{u, 0) = m\R{u, 0) ^ 0. 


Hence, 


h'mW = 


d(^Q{uj,0) ml 


dF{u)-,B) 

dB 


{V) 


= -K{u,B;V). 


(4.5) 


R{uj,0) df^F{uj,B) 

Denote by P{zX) ■= {z — u)^ + hi{(){z — u)^~^ -\ - \-hm{C) the Weierstrass polynomial. Since 

R{z, () 7 ^ 0, the zeros of Q and R in x ^^^(O) coincide. It is well known that for small enough 

S G (0, (5i] and |C| < 5, there are exactly m 2 ;-roots (counting multiplicity) of P{zX) = 0 and these 
roots are given by a one or more convergent (and possibly multivalued) Puiseux series. We denote 
these m branches of (one or several) Puiseux series by ^j{C). 

The method of the Newton diagram (see e.g. [44] and [38l Theorem XII.2]) implies that, if the 
roots Dj(C) are given by more than one Puiseux series, then h(^(0) = K{oj, B,V) = 0 and all the 
roots have the asymptotics ^j{C) = w + o(C^'^™'). 

In the case when the roo ts D ,(C) are given by one P uise ux series, this series can be written in 
powers of and formula (4.4) follows from equalities (4.5) and (Dj — a;)™' -|- + 

■ ■ ■ + hm{C) = 0) see details in the proof of [2T1 Proposition 3.5]. □ 


5 Variational characterization of the nonlinear spectrum 


Recall that is the indicator function of 
solution to 


A function y G [ 01 ^ 02 ] is called a nontrivial 


y'fx) = -uj‘^y{x) hi{x) + [h 2 {x) - hi{x)]xrAy'^{x)) 


(5.1) 


if y ^ 0 and (5.1) is satished for a.a. x G [ 01 , 02 ]. If problem ( |5.1[ ), (1.2), (1.3) admits a nontrivial 
solution y, we say that y is an eigenfunction and uj is an eigenvalue of (5.1), (1.2), (1.3). In short, 
we will say that u is a nonlinear eigenvalue. The set of nonlinear eigenvalues is denoted by 
Note that 

is symmetric w.r.t. iM. (5.2) 

Indeed, if co = iC G and y is an associated eigenfunction, then for B = [bi + (62 — 5i)yj,^ (|/^)] one 

has u G E(R). Since coi = iC G E(i?) and y is an eigenfunction associated with uji (see Lemma 3.3 
(ii) an its proof), we conclude that uji is a nonlinear eigenvalue with yi '.= iy as an eigenfunction. 
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Let us notice also that B = \bi + (62 — belongs to A and, moreover, belongs to the 

set ext A of extreme points of A. So, 


C S[extA] C S[A]. 


(5.3) 


Let K be the functional defined in Proposition 4.4| (ii) and let cone iL(a;, S; A 
nonnegative convex cone generated by the set K{u, B; A — B) = {K{u, B;B — B) : 


— B) he the 
B e A}. 


Theorem 5.1. Let u ^ iM. Then the following statements are equivalent: 

(i) u is nonlinear eigenvalue, 

(ii) there exists Bq & A such that oj & T^Bq) and cone K{u, Bq; A — Bq) ^ C. 

such that 


If (i)-(ii) hold true, then there exists an eigenfunction y of (5.1), (l.lS), (1.5 


Bo(x) ^ + [62(x) - 6i(j:)]Xj.^(y^(x)) a.e. on {a,.a2). 


(5.4) 


The proof (see Section 5.2) relies on the rotation properties of 6 described in Section 5.1 


5.1 The turning interval and rotational properties of modes. 

By nrgz we understood the multivalued argument function of a complex variable 2 ; 7 ^ 0. In this 
subsection, we write 9{x) for 9{x,u;B). 

Lemma 5.2. Let B E L^ and u G T,{B). Assume that a = Reca > 0 and B{x) > 0 a.e. on ( 01 , 02 ). 
Then: 

(i) There exists a unique subinterval of [ 01 , 02 ] with the properties that 


9{x)dx9{x) G C_ when x < x^,, 9{x)dx9{x) G M when x G [x^:,x*], 

and 9{x)dx9{x) G C+ when x > x*. (5.5) 


(It is supposed that x* < x* and that, in the case x* = x* , the turning interval [x*,x*] degenerates 
into a single turning point x* = x*. The meaning of ’turning’ is explained by statement (5.1) below). 

(ii) If x^ < X*, then 


B = 0 a.e. on (x*,x*) and dx9{x) is constant on [x*,x*]. (5.6) 

In particular, B > h a.e. yields x* = x*. 

(Hi) There exists at most one point xq G [ 01 , 02 ] such that 9 {xq) = 0. If xq does exist, it belongs to 
[x*,x*]. (In the sequel, when Xq does not exist, we assume {xq} = 0.) 

(iv) The multifunction arg6'(x) has a branch arg^ 6*(x) that is defined on [ 01 , 02 ] \ {xq} and has the 
following properties: 


arg^ 6 *(x) is differentiable on [01,02] \ {xq}, 

dx arg^ 9 < 0 if X < x^, dx arg^ >0 if x > x*, dx arg^ 0 = 0 if x E [x*, x*] \ {xq}. (5.7) 

(v) dx arg^ 0 (x) is bounded on [oi, 02 ] \ {xq}. In particular, for every G [—vr, vr), 

the set {x G [01,02] \ [x*,x*] : arg 0 (x) = .^(mod 27 r)} is at most finite. 
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Proof, (i)-(ii) Since Imo; < 0 and Reca > 0, we have Imcj^ < 0 and 

dx^^{0dx9) = \m.{9d^9) = \xn.{—9upB9) = B\9'^ \m{—uP‘) > 0 a.e. 


(5.8) 


For X = ai and x = 02 using (3.1) and (1.3), one gets 


Im[6'(ai)9a;6'(ai)] = ui Ini(—ica) = —uia < 0 


and 

U 2 — ica 

Hence, lm{9dx9) is a nondecreasing continuous function having a zero in [ 01 , 02 ]. Denote 


X* = niin{x : \m[9{x)dx9{x)\ = 0} and x* = niaxjx : \m[9{x)dx9{x)\ = 0}. 


\m[9{a2)dx9{a2)] = Im 


— IOJV2 


Then 


lm[9{x)dx9{x)] < 0 for X < X*, and \m[9{x)dx9{xy\ > 0 for x > x*. (5.9) 

3S 

Ini[ 6 *(x)cl 3 , 6 *(x)] = 0 for all x € [x*,x*]. (5.10) 


When X* < x*, (5.8) implies 


Combining this with 9^0 = —u^B9, we get (5.6). 

(iii) It follows from (i) that every x-zero of 9 belongs to [x*, x*]. If there exist two distinct zeroes, 
then ( 5.6| ) yields dx9 = 0 on [x*,x*] and so 0 = 0 on [ 01 , 02 ], a contradiction (to (3.1)). 

(iv) Consider an arbitrary continuous on [ 01 , 02 ] \ {xq} branch arg^ 6 '(x) of arg 6 '(x). Combining 
the equality 


lm{9dx9) j 9x9 

— 1 ^ 1 ^— = Im — = Im da, In 9 = dx arg^ 9 for x / xq. 


(5.11) 


with ( 5^ and (5.10) we get (iv). 

(v) Due to formula (5.11), it is enough to prove that dxa.Tg^, 9 is bounded in a punctured neigh¬ 
borhood of xo assuming that xq exists. 

To be specihc, consider the case x > xq. Since 9{xo) = 0, one has 9'{xo) 7 ^ 0 and 9{x) = 
9'{xo){x — Xo) — Ixo^^ ~ s)9{s)B{s)ds. Hence, 


Im(6'da;6') 

dxaTgJ= = |o;| 



Xo ^Xo 


= \0J\ 


, , , , (x — s)9(s)9{t) 

B{s)B{t) Im--ds dt = 



XQ JXQ 


^ Bis)B{t) Im d. *. (5.12) 

\9{xW 


where p{t,s) lies between t and s. Since d jg bounded in the square (s,f) G (xo,x)^ (recall 

that 9'(xo) 7 ^ 0) and 9,9' G C'[oi, 02 ], we see that the last integral in (5.12) is uniformly bounded for 
X > Xq. □ 

Now consider modes 9 corresponding quasi-eigenvalues on the axis iM. In this case, the analysis 
is slightly different, but simpler, since 

9(x, z; B) are real when 2 ; G iM and B G ^^(ai, 02 )- (5.13) 

Lemma 5.3. Let u G D(i?), u = i(3, /I G M, and B{x) >0 a.e. on (oi, 02 ). Then: 
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(i) There exists a unique subinterval of [ 01 , 02 ] such that 


6{x)dx6{x) <0 if X < x^, 


9{x)dx0{x) = 0 if X E [x^^x*], 

and 9{x)dx9{x) >0 if x > x*. 


(a) If x^ < X* , then i?(x) = 0 a.e. on and 9{x) is a nonzero constant function on [x^^x*]. In 

particular, if B{x) >0 a.e. on ( 01 , 02 ), then x^ = x*. 

(Hi) If 9{xo) = 0, then xq = x^ = x*. 


Proof. Let x^ = min{a; G [ 01 , 02 ] : 9{x)dx9{x) = 0} and x* = max{x G [ 01 , 02 ] 
The existence of x* and x* follows from (1.2)-(1.3) and 9dx9 G C[oi, 02 ]. Since 


9{x)dx9{x) = 0 }. 


dx{9dx9) = {dx9Y + 9819 = [8x9)“^ + I3‘^B9‘^ > 0 a.e. on (oi, 02 ), 


we obtain (i). The eqnality 

px* px* 

0 = 9{x*)8x9{x*) — 9{x^:)8x9{x^:) = / 8x{98x9) dx = / [{8x9)“^ + (3‘^B9‘^]dx 

J X:^ J X* 

implies 8x9{x) = 0 and B{x) = 0 a.e. on [x*,x*]. When x* < x*, this and 9^0 yield statement (ii). 
Statement (hi) follows from (ii). □ 


5.2 Proof of Theorem 15.11 


Let a = Re CO 7 ^ 0. Due to Lemma 3.3 (ii) and (5.2), without loss of generality one can assume a > 0. 
Suppose u G T‘{Bq) for certain Bq G A. Put 

8x9{a2, uj, Bq) 


Cq = 




By (4.3), 


a 


dF{u.; B,) ^V) ^ 


Define the set So C C by 

8 F{u-Bq) 


Sq := Cn 


8B 


8B 


[A - Bq] = 


r(s,co;Ro) R(s) ds 


fai 


^a2 


9\s,u;Bq) R(s) ds : R G 


Br 


'ai 


The definition of K{uj, Bq] V) in (4.4) implies the equivalence 

coneRoT^C cone Kipj, Bq] K —Bq) 


(5.14) 


(5.15) 


Proof of implication (i) ^(ii). Let y be an eigenfunction of (5.1), (1.2), (1.3) associated with 
CO. Let us dehne Bq by equality (5.4) for all x G [ 01 , 02 ]. Then co G H^Bq). More over , there exist 
C G M \ {0} and G [— tt, vr) such that y{x) = Ce~^^*^‘^9{x,u] Bq). It follows from (5.4) that 


DC 'i _ / ^ 1(^)1 when 9^(x) G e*^*C_ 

° ^ \ ^ 2 ( 2 :), when 6 *^(x) G e*^*C+ 


(5.16) 


This implies that for every V G A — Rq the following assertion hold: 

R(x) > 0 for a.a. x such that 9‘^{x) G e*^*C_, 

V{x) <0 for a.a. x such that 9‘^{x) G e*^*C+. 
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So 9‘^V ds G e*^*C_ for each 1/ G A — Bq. Thus, cone So C e*^*C_. Equivalence (5.15) concludes 
the proof of (ii). 

Proof of implication (ii) ==^(i). Put arg^ 6 '^(a;) := 2arg^6'(a;). Define the sets Ef C [ 01 , 02 ] by 

Et ■= suppffea - Bq] + , E- := supp[5o - 5i]+, 

(for [■]+ see the notation part of Section [^. 

Lemma 5.4. Let Reco > 0 and 9{x) = 9{x,u; Bq) . Then 

9^[E+] U {-9^[E;]) C cone5o, 

i.e., the image of the set Ef under the function 9"^ and the image of the set E~ under the function 
{—9"^) are subsets of the closed cone cone So 


Proof. Recall that the point xq was defined in Lemma |5.2| and that, in the case when xq does not 
exist, the notation {xq} means the empty set. 

Step 1. If Xq exists, 9{xq) = 0 G coneRo- 

Step 2. Suppose Xi G Ef \ {xq}. Then Vf^{x) := X(xT_-e,xi+e){,x) [62 — Bq]^{x) has a nonempty 
support for every £ > 0. Let us show that (for a suitable choice of values of arg) 


arg 


■ dE{u-BQ) . 

dB > 


arg^ 9‘^{xi) —)• 0 as e —>■ 0 . 


By Lemma 5.2 (iv). 


arg^ 9‘^ is dehned and continuous in a vicinity of xi. 


(5.17) 


(5.18) 


For small enough e > 0, | arg^ 9'^{x) — arg^ 9'^{xi)\ < 11/2 in the interval x G (xi — £, xi + e) and 
therefore the equality 


dE /■“■ 

- i 


02 ds = 


I'Xi+e 


9^[b2-BQUds, 


' Xi—£ 


implies that CQ^{Vf~) G cone02[(xi — £,xi + e)] \ {0}, where 6*^[(xi — e,Xi + e:)] is the image of 
(xi —e,Xi Re) under 0^. This and (5.18) yields (5.17). 

Step 3. Suppose xi G E~ \ {xq}. Then, similar to Step 2, Vf{x) := —X{xi-e,xi+e){.x)[BQ — 6i]+(x) 
has a nonempty support for small e and arg [iPoff (K~)] ~ ^^(xi) — tt —)■ 0 as £ —>■ 0. In other 

words, 

dE 

argCo^(K) - arg[-0^(xi)]0 as e0. 

Step 4- Since B + Vf" G A, we see from Steps 2-3 that 6*^(xi) G cone Rq whenever xi G Ef \ {xq}, 
and [—6'^(xi)] G coneRo whenever xi G E~ \ {xq}. This completes the proof. □ 

Lemma 5.5. Let Rea; > 0. Assume that coneRo A contained in a certain closed half-plane 
where G [—7r,7r). Then: 

(i) eguality holds for a.a. x G [ 01 , 02 ], 

(ii) a; G and y = is an associated eigenfunction of (5.1), (1.2), (1.3). 
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we see 


Proof, (i) Recall that the set E is dehned by (2.4). From the dehnition of E and Bq G 
that E C Xq U Ef U E~ , where Xq is of zero measure. Clearly, cone So C e*^*C_. By Lemma 5.4 

e‘^[Et] C e*«‘Cr and C 

So Ef n Ej C {x : 9‘^{x) e By Lemma 5.2 (iv)-(v), 

the set Xi := {x : 9{x) G \ [x*,x*] is at most hnite. 


(5.19) 


(5.20) 


Thus, E <Z XqU XiU [x*,x*] U {Ef \ E~) U {E~ \ Ef). On the other side, Bo{x) = bi{x) = 0 a.e. on 
[x*,x*] due to Lemma 5.2 (ii). So X 2 := ([x*,x*] O E)\ {Ef \ E~) is of zero measure. 

Summarizing, we see that the interval [ 01 , 02 ] can be written as 

[oi, 02 ] = ([oi, 02 ] \ E) U Xo U Xi U X 2 U {Ef \ Ef) U {Ef \ Ef). 

On each of the sets of positive measure in the right side, equality ( 5.16| ) is fulhlled for a.a. x. 

Additional explanations are needed for fact that 9'^{x) G e*^*C+ for a.a. x G E~ \ Ef. Assume 
that xi belongs to the se t X 3 consisting of x such that x G E~ \ Ef and 9‘^{x) ^ e*^*C+. By (5.19), 
9^{xi) G e*^*]R. By (5.20), X 3 C Xi U [x*,x*]. Assume that measX 3 > 0. Then meas(X 3 n [x*,x*]) > 
0. This yields that meas [supp[i? — 61 ]+ O (x*, x*)] >0. The latter contradicts Lemma 5.2 (ii). This 
contradiction shows that measXs = 0 . 

(ii) Equality (5.16) implies (5.4). The latter yields (ii). □ 


Equivalence (5.15) and Lemmata 5.4, 5.5 imply statement (i) of Theorem 5.1 


The last assertion of Theorem 5.1 follows from Lemma 5.5 (ii) and its proof. 


6 Local extremizers and optimization for a fixed frequency 


6.1 Definitions of various optimizers 


We say that a G M is an admissible frequency if a = Reca for some admissible quasi-eigenvalue u. 
So ReS[A] is the set of admissible frequencies. Its properties are considered in Subsection 6.2 and 
Appendix |A.3[ where it is proved, in particular, that high enough frequencies are admissible at least 
for most popular settings of quasi-eigenvalue optimization problem. 


Definition 6.1 ([23])- Let a be an admissible frequency. 

(i) The minimal decay rate d min fo) for the frequency a is defined by 


/dmin(«) := inf{/? G M : a — i/d G S[A]}. 


(ii) If ca = a — id min (o) is a quasi-eigenvalue for a certain admissible structure R G A (i.e., the 
minimum is achieved), we say that u and B are of minimal decay for the frequency a. 

A pair (cn, R) G C x is called admissible if R G A and u G S(R). 

Definition 6.2. We say that an admissible pair (cuo, Ro) is a local maximizer oflmoj for the frequency 
Recuo if there exist £ > 0 such that Imcui < Imcuo for any admissible pair (a;i,Ri) satisfying the 
following three conditions 


Retui = Recuo, \oJi — cuol < £, 


-Bi — RqIIi < 
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P 


2[A]nP = AjUA2UA3UA^ 


Re5)[A] 3 ReAj^U ReA 2 U Re Ag U ReA^ 


Figure 1: A schematic drawing (see explanations in Remark 6.2) of the part of the set of admissible 
quasi-eigenvalues S[A] in the rectangle T) = {Re 2 ; G (0,1-2), Im^; G (—0.015,0)} and of intervals 
of admissible frequencies generated by E[A] for the following choice of parameters: z/i = 1, 


z/2 = -l-C)0, hi = 90, 62 = 110, 02 — Oi = 1. 

mark one of local minimizer for /^min and the corresponding quasi-eigenvalue; 
’x’ mark quasi-eigenvalues of minimal decay for frequencies a in ReA 2 ; 

’o’ mark co-parts of local minimizers of Imco for frequencies a in Re A 3 . 


Local minimizers oflmu for a particular frequency are dehned in a similar way. A pair (co, B) is 
a local extremizer if it is a local minimizer or a local maximizer. 

An admissible frequency is a local minimizer for /?min if there exists e > 0 such that d m^n fan) < 
/3min(tt) for all admissible a in {ao—e, aQ+e). If for certain e > 0 and all admissible a from a punctured 
neighborhood (ccq — £, ao) U (ctoi cto + ^) fhe strict inequality d min fan) < /5min(tt) holds, ccq is said to 
be a strict local minimizer. 

Usually, in applied and numerical literature, local minimizers of Imco are considered without 
frequency restrictions. It seems that such approach goes back to the dehnition of [H]. Below we give 
a rigorous adaptation of that dehnition suitable for the optimization over A. 

Let us dehne on the set of admissible pairs (cn, B) the decay rate functional Dr(ci;, B) := — Imcn. 

Definition 6.3. An admissible pair (cuo, Bq) is said to be a local minimizer for Dr if there exist 
£ > 0 such that Dr(a;o,i?o) < E)r(a;i,i?i) for every admissible pair {oji^Bi) with |a;i — cuol < £ and 
ll-Bi — -Bolli < £. Local maximizers for Dr are dehned similarly. 


Remark 6.1. Clearly, if {ooq^Bq) is a local minimizer (maximizer) for Dr, then {ujq^Bo) is a local 
maximizer (resp., minimizer) of Imo; for the frequency Recuo- 


Remark 6.2. Figure is drawn to illustrate the above dehnitions of various optimizers. It also 
provide a rough impression about a part of S[A] in B for the one-side open case z/i = 1, z /2 = -|-oo 
(cf. Example 6.2). The drawing is partially based on the MATLAB computations of the nonlinear 
spectrum (see Section 11 and Figure]^ and partially on Proposition 6 . 6 , Theorems 


8.1, 8.2, 9.1 


and formula ( 10 . 1 ). 

In agreement with Theorem |8.2| E[A] PV consists of four ‘clouds’ Ai, ..., A 4 . These ‘clouds’ are 
connected components of E[A] produced by all possible perturbations (inside A) of quasi-eigenvalues 
(Uo,! = i^o(^ 2)5 1 ^ 0,2 = ..., Wo ,4 = ^ 3 ( 62 ), resp., which are generated by the homogeneous 

structure Rq = ^ 2 , see (3.7[) and Section]^ 
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The projections of these components to the real line ReAj are snbsets of the set of admissible 
freqnencies ReS[A]. However, we do not know whether [J^^^ReAj are all admissible frequencies in 
the interval (0,1.2). Indeed, it is not known if there exist any admissible quasi-eigenvalnes u with 
Reca G (0,1.2) lying ontside of V. One of the difficnlties in the stndy of a; G S[A] with large |a;| is 


that, globally in C, H(i?) is not a continnous set-valued function of B, see Remark 4.1 


In these setting, the frequency 0 is not admissible, see Remark 7.2 


6.2 Admissible frequencies and examples of admissible families 


Recall that the notation [a:] ([a:J) stands for the ceiling (floor) function. 

In this subsection we consider several most reasonable and popular admissible families for which 
the choices of the constraints has a Physics motivation. For all this families we provide a quantitative 
version of the statement that high enough frequencies are admissible. These estimates on the ranges of 
admissible frequencies are quite technical. They are derived in Appendix |A.3| together with the proof 
that high enough frequencies are admissible in more general situation, which is described Proposition 
6.5 (iii). This proposition also considers other properties of the set ReS[A] of admissible frequencies 


and, roughly speaking, shows that the set ReS[A] cannot have some ’wild’ structure. 

Example 6.1. Let the constraint functions 61^2 be constant and 0 < 61 < 62 - Assume that ui = i >2 =: 
u > 0 and equals either 61 , or 62 - These settings represent the admissible family associated with 
two side open optical cavities of length 02 — ai consisting of transparent dielectrics with permittiv¬ 
ities in the range from bi to 62 |12]- At the endpoints 01 ^ 2 , the cavity is connected to half-inhnite 
homogeneous outer media consisting of the dielectric material with permittivity In Optical En¬ 
gineering modeling, the outer medium is usually the same as one of dielectrics with extreme allowed 
permittivities (usually, either 61 = = 1 represents vacuum, or 62 = corresponds to a cavity 


connected to a waveguide of permittivity 62 0)- Under this choice of parameters, formulae (|A.9[) 


and (A.13) of Appendix A.3 show that 


all frequencies in the range |a| > 


TT 


(02 - ai)b 


1/2 


K 


1/2 


l 1/2 ,1/2 




are admissible. 


Example 6 . 2 . Let the constraints 61^2 be constants and 0 < 61 < 62 - Assume that uf = bi and z /2 = 00 . 
One of Optics interpretations of these settings is one-side open cavity, which, in particular, is used 
for the laser modeling |13]. The Dirichlet boundary condition 1 /( 02 ) = 0 generated by the parameter 
z /2 = 00 corresponds to a perfect mirror (a plate of a perfect electric conductor perpendicular to the 
x-axis). So EM waves are reflected without dissipation from the interface plane x = 02 . Through the 
interface plane x = ai the waves radiate into homogeneous half-inhnite medium with permittivity 
uf. It was proved in [22] that 


Example 6.3. In Mechanics models involving the equation for transverse oscillations of a nonhomo- 
geneous string, the linear density B{x) of the string on some intervals may be much less than on 
others. In such cases, the massless string approximation is employed. This means that B{x) (and 
so the lower constraint bi{x)) can be equal to 0 on certain intervals. Keeping this in mind, consider 
constant constraints 61^2 £ [0, -l-cxo) satisfying 0 < 61 < 62 - Assume that z/i = 0 and 0 < 1/2 < -|-cxd. 
The choice z/i = 0 leads to the Neumann boundary condition y'{ai) = 0 and corresponds to the 
assumption that the left end of the string is ’’free” (i.e., a massless ring at the end Oi of the string is 
sliding without friction on a pole pointing in the transverse direction). The Mechanics interpretation 


all frequencies in the range |a| > 


TT 


(02 - ai)b 


1/2 \ 2 


1 

X + 


K 


1/2 


l 1/2 ,1/2 


K' 


^ are admissible. ( 6 . 1 ) 
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of the condition (1.3) at 02 is the linear damping condition (which arises, for instance, if a massless 
plate immersed in a flnid is connected to the right end 02 ). In snch settings, z /2 is the damping 
coefficient (if the units are chosen such that the tension of the string equals 1). In this context, V 2 
has no natural connection with extreme allowed densities of the string. The estimates of Appendix 
lA. 31 shows that for this model 


all frequencies in the range |q!| > 


37r 


2(02 - ai){hy‘^ - h\^‘^) 


are admissible. 


Another interpretation of the Neumann condition y\ai) = 0, which makes also sense from Optics 
point of view, is given in the next example. 


Example 6.4. Let us assume additionally to the settings of Example |6.1| that 02 = —Oi and let us 
restrict optimization to structures symmetric w.r.t. origin, i.e., to B such that B{x) = B{—x) and 
hi < B <1)2 a.e.. This model of a symmetric resonator is used often to reduce computational costs, 
see e.g. [19]. Since B is an even function and vi = z/ 2 , it is easy to see that modes y of such a 
resonator are either even, or odd, and so either condition |/'( 0 ) = 0 , or the condition |/( 0 ) = 0 is 
fulhlled for them. Thus, the quasi-eigenvalue problem and the related optimization problem can be 


essentially reduced to two: the problem considered in Example 6.3 and a problem similar to that of 


Example 6.2 (possibly, with different vi). 

Denote 6 ™®'^ := esssup 6 i(a;) and := ess inf 62 (a:)• Recall that intervals are, by dehnition, 
connected subsets of M and that an interval called degenerate if it consists of one point. 


Proposition 6.5. {i) If I is a connected component o/ReE[A], then either I is a nondegenerate 
interval, or I = {0}. In particular, the only possible isolated admissible freguency is 0. 

(ii) Suppose fc™*™ < 6 ™™. If b'f^^ < vi or > z/ 2 , then 0 is an admissible freguency. 

(Hi) Suppose 6™®'^ < Then high enough freguencies are admissible. Moreover, 


Re E[A] D < |a| > 


27r 


(o2-oi)(y6=-y6r 


( 6 , 2 ) 


Proof, (i) It is enough to show that if lvq G T,{Bq), Bq G A, and Reozo 7 ^ 0, then Bq can be piecewise 
linearly perturbed in A such that perturbations of uq do not preserve Reo/o- 

Taking nontrivial V G A — Rq consider perturbations Bo + (V with sufficiently small C > 0. If for 
perturbed quasi-eigenvalues u{() the real part is preserved, we hx a certain branch and certain 
sufficiently small (^1 G (0,1) such that Bi = Bq + (iV is not an extreme point of A. By Theorem 
5.1 and the remarks before it, cone A'(a;(Ci), Ri; A — Ri) = C. By Proposition 4.4, a;(Ci) can be 


perturbed with change of the real part (see also Proposition 10.2). 


(ii) follows from Proposition 3.4 (ii). 


(iii) For the settings of Example 6.2 
in slightly different form. 


( 6 . 2 ) follows easily from formula ( 6 . 1 ), which was proved 


m m siigntiy amerent torm. The proof in general settings can be obtained using the same 
idea, but requires more technical considerations and involves a number of different cases of mutual 
arrangements of z/i_ 2 , 


and 6 ™™. This proof is given in Appendix A.3 


□ 


6.3 Existence of optimizers and properties of E[A] 

Proposition 6.6. The set S[A] is closed. 
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Proof. Since A C the proof requires some modifications of standard compactness arguments. For 
each i? G A there exists unique H{x) = H{x] B) such that 

B{x) = {b 2 {x) — hi{x))H{x) + 6 i(x) a.e. on (oi, 02 ) 

and H{x) belongs to the set A* of L“-functions satisfying 

0 < H{x) < 1 a.e. on (ai, 02 ) and H{x) = 0 a.e. on (oi, 02 ) \ E. 

Let us take R > ||& 2 ||i- One can see that 


the map iL(-; B) 1 —>■ B{-) is a bijection of A* to A continuous 

from the weak* topology of to the metric topology of 


,p) (6.3) 


(see Subsection 4.1, the continuity is sequential and w.r.t. the corresponding induced topologies on 
A* and A). Lemma 4.1 (i) and the sequential Banach-Alaoglu theorem applied to A* complete the 
proof. □ 


Proposition 6.6 and S[A] C C_ imply three following results. 


Corollary 6.7. (i) For each admissible frequency a, there exist a structure B and the quasi¬ 
eigenvalue oj of minimal decay, i.e., u = a — i/5niin(a) belongs to S(i?) for certain B E A. 

(a) The pair {pj, B) from statement (i) is a local maximizer oflmu for the frequency a. 


Proposition 6 . 8 . The function is lower semicontinuous, i.e., i/{an}(f C ReE[A] and an —)■ ao, 
then /3min(ao) < liminf/3min(a„). 


Proposition 6.9. If an admissible frequency a is a local minimizer for (drain, then the pair {pj, B) 
from Corollary 6.1 is a local minimizer for Dr. 


The existence of local minimizers for (drain and Dr is considered in next two sections. 


7 Non-uniqueness and local minima for small a 


Pure imaginary quasi-eigenvalues ca G iM are the subject of special interest in Mechanics, see e.g. 
[SI El- They corresponds to overdamped and critically damped oscillations. 

Proposition 7.1. IfReoJo = 0 and {uo,Bq) is a local extremizer oflmoj for the frequency 0, then 
either Bq = bi, or Bq = 62 - 7n particular, if Bq is a structure of minimal decay for the frequency 0 , 
then either Bq = hi, or Bq = 62 . 


Proof. Assume that (wq, Bq) is a local extremizer of Imca for frequency 0 and that Bq ^ bi, Bq ^ 62 
(in the sense of L^-space). Then it follows from (4.3) and Lemma 5.3, that there exist E A — Bq 
such that > 0. This, Proposition 4.4, and the symmetry of E(R) w.r.t. iM imply 


dFpjQppp I 

dB 

that there exist V-f G 


— Bq and a.;±(C) G E(i?o + C^ 2 ^) such that for small enough C > 0, 


Reci;±(C) = 0, ±[Ima;±(C) — Imcuo] > 0, and Ci;±(C) —)■ cuq as C —)■ 0 

(the arguments are similar to [211 Section 4.4], note that and Vff are not necessarily different). 
Thus, (a;o,i?o) is not a local extremizer. □ 

Theorem 7.2. Assume that a = 0 is an admissible frequency. Then 
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(i) Sim(6i, 62 ) := ilR n [S(6i) U S(& 2 )] is a nonempty closed set. 

(ii) /3min(0) = min |Ima;| 

weSini(6l,fe2) 

(Hi) a = 0 is a strict local minimizer for /?min whenever co = —i/3mm(0) does not belong to Smuit(&i) U 
Smuit(^ 2 ) (i-e., whenever the quasi-eigenvalue of minimal decay u is not a multiple quasi-eigenvalue 
of bi and 0 / 62 ^- 


Proof, (i) and (ii) follow from the assumption 0 G ReE[A], Proposition 7.1, and Proposition 6.6 
If 0 is an isolated point of ReE[A], statement (iii) follows from the dehnition of a strict local 
minimizer for (see Section 6.1). Let us denote by I the connected component of ReE[A] 
containing 0, and let us prove (iii) in the case when I 7 ^ {0}. 

(i), / contains a nonempty symmetric interval (—a!o,Q!o). 
Then also = lim infQ,_>o+ /dn 


By Lemma 3.3 


11 


and Proposition 6.5 
Consider /?* := lim infQ,_^o-A 


mmvCtJ- J-iien also p* = iimimQ,_>o+Pminl«) due to Lemma 3.3 (ii). 
If A* = + 00 , it is obvious that a = 0 is a strict minimizer. Consider the case A* < + 00 . 
Since E[A] is closed, there exist sequences C A and oj±n G E(i?„) such that 


\m.uj±n —>■ —A*) Re Win —)■ 0 ± 0 , and 

the sequence Bn converges to a certain R* G A w.r.t. the metric p of Subsection 4.1| 


(this can be proved with the use of the map H{--,B) 1 —)■ R(-) from (6.3) and the sequential Banach- 
Alaoglu theorem applied to A*). 

Let us s how that a;* := —iA* is a multiple quasi-eigenvalue of R*. Indeed, ca* G E(R*) due to 
Lemma 4.1 (i). One can take £ > 0 such that De(a;*) does not contain other points of E(R*). For 
large enough n, the disc ©(.(ca*) contains at least two distinct quasi-eigenvalues u±n of R„. Hence, it 
follows from Proposition |4.2 that the multiplicity of ca* as a quasi-eigenvalue of R* is > 2. 

Now, the assumption that u = — iAmin(O) is a simple quasi-eigenvalue of any corresponding 
structure R of minimal decay (by Proposition 7.1, R is either bi, or 62 ) yields 00 ^ 00 (since they 
have different multiplicities). So A* > Amin(O). This completes the proofs of statement (iii). □ 


Now we will show that each of the statements (i)-(iii) of Theorem 7.2 takes place for constant 
constraints 61^2 under the assumption that at least one of the constants ^ 2 ^^ is not in the interval 

Recall that the positive function Ki{h) dehned in Proposition |3.4| for positive b is given by Ki(b) = 
Ki{b, ui, U2) : = 


l+l^l/^2 


''1 

61/2’ 


G /2 

"2 


Dehne 


K 2 {h) = K 2 {h,vi,V 2 ) ■=s.gn[l-Ki{h)] 


K,{h) 


i/Vh 


for 6 > 0 , 


l + iFi( 6 ) 

772 ( 0 ) = 772 ( 0 , z/i, 1 / 2 ) := exp(—2/z/i — 2 /Z/ 2 ) (when z/i = 0, we suppose 772(0) := 0). 


Corollary 7.3. Let hi and &2 be nonnegative constants such that b^^ < ui or > z/ 2 . Then: 

(i) a = 0 is a strict local minimizer for Amin and Amin(O) = ~ 2(02-01) i^3'X{772(5i), 772(62)}- 

(ii) 7 / 772 ( 61 ) = 772 ( 62 ), then the set of structures of minimal decay for the frequency 0 consists of two 
structures, 61 and 62 . 

(iii) 7 / 772 ( 61 ) > 772 ( 62 ), then hi is the unique structure of minimal decay for the frequency 0 . 

(iv) 7 / 772 ( 61 ) < 772 ( 62 ), then 62 is the unique structure of minimal decay for the frequency 0. 


Remark 7.1. By Proposition 3.4, the quasi-eigenvalue of minimal decay Uq = — iAmin(O) has mul¬ 
tiplicity 1 for every associated structure of minimal decay in all the cases of Corollary 7.3 (cf. 
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Proof of Corollary \7.3\ Proposition |3.4| implies that 0 is an admissible frequency. It follows from 
(|3.8[) that K 2 {b) < 0 exactly when S( 6 ) fl iM = 0 . On the other side, 7 ^ 2 (&) > 0 exactly when 


S( 6 ) n iM = {a;o(&)}, where uJo{b) is dehned by (3.7). In this case, K 2 {b) = exp ( 2 (a 2 — Oi) Ima;o(&))- 
This, Theorem 7.2, and Proposition 3.4 complete the proof (see also Remark 7.1). □ 


Example 7.4 (non-uniqueness of optimizer). Let C G M+. Put bi{x) = b 2 {x) = 40^, z/i = z /2 = 

Oa/ 3. This correspond to optimization problem for a cavity constrained by < B(x) < 40^ and 
placed into the outer medium (—cx),ai) U (a 2 ,-|-cxD) with permittivity (cf. Example 6.1). The 
length 02 — Oi is not important for the present example. 


Let ujq := 


2C{a2—ai) 


' In (7 -|- 4 a/ 3) . Then the pairs (cuo, &i) and (coq, ^ 2 ) are: 


(a) local minimizers for Dr and local maximizers of Imo; for the frequency 0. 
Moreover, 

(b) 61 and 62 are structures of minimal decay for the frequency 0 , 

(c) the frequency ao = 0 is a strict local minimizer for /?min- 


These statements follow from Corollary |7.3| (i)-(ii) and straightforward calculations. This example 
answers partially the question of uniqueness of optimal structure discussed in [HI |23l |20] 

Remark 7.2 (cf. HE]). Considering the case when constant functions 61^2 do not satisfy the assump¬ 
tion of Corollary 7.3, we obtain the following result: if vi < B[x) < V 2 a.e., then S(i?) fl iM = 0. 

imply that 0 ^ ReS[A] when vi < b\^'^ < b^^ < V 2 - 


Indeed, Theorem 


7.2 


and Proposition 


3.4 


8 Local minima in the case of small contrast 


Local minimizers for /dmin and for the decay rate functional Dr were dehned in Section 6.1 


Theorem 8.1. Let b G M+ \ [uf, z/|]. Then for every natural number N there exists e > 0 such that 
for any admissible family A defined by constant constraints 


bi = b — ei, b 2 = b + 82 with t) < si -\- 82 < s and £ 1^2 > 0 , 


the following statements hold: 

(i) In the interval |tt| < exist at least 2N -|- 1 local minimizers for /?min- 

(a) There exist at least 2N + 1 local minimizers and at least 2N -|- 1 local maximizers for Dr with 
associated frequencies in the interval |a;| < 


Remark 8.1. In the case b G (z/^,z/|), 
extremizers to 2N everywhere and the interval to |a| < 


TVtt 


Theorem 8.1 becomes valid if one change the number of 

Figure j^i 


bl/2(a2—ai) 


see 


and Remark 


6 . 2 ). 


The theorem follows from the next (more general) result, its proof, and Proposition 3.4 
Recall that a set S' C C is called path-connected if for any two points zi ^2 ^ S' there exist a 
continuous function / : [0,1] —>■ S such that /(O) = Zi and /(I) = Z 2 - 


Theorem 8.2. Let Bq G L^ , Bq{x) > 0 a.e., and let T> he an open hounded subset of C. Assume 
that E(Ro) n ^ consists of n E N distinct points {a;oj}"=i and there are no points of E(Ro) on the 
boundary hdV. 
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1) Then there exists e > 0 such that for any functions &i,2(a^) satisfying 

0 < 6i < i?o < &2 o-e. and 0 < ||6i — 62II1 < ^ (8.1) 


the following statements hold: 

(i) S[A] n V consists of n (disjoint) closed path-connected components Ai, ..., A„ satisfying 
uoj e Aj, j = 

(a) There exist at least n local minimizers and at least n local maximizers 

Bf^^) for Dr such that ^ Aj, 3 = 1 , ..., n. 

(Hi) For each a G U^ReAj, there exists at least one local minimizer and at least one local 
maximizer oflmu for the freguency a. 

2) If, additionally, V has the form {ai < Rez < 02, —(3o < Imz < 0} with certain ai,2,/9o ^ 

and if e is chosen as in statement ( 1 ), then in the freguency interval {ai,a 2 ) there exists at 
least one local minimizer uq for (drain- This minimizer Oq can he chosen such that /3min(«o) = 
min^6S[yi]n© | Imw|. 


The rest of this section is devoted to the proof of Theorem 8.2 


Proposition 8.3. Assume that i?i,i?2 G and a bounded open set D C C are such that for all 

te [0,1], 

'E{B{t)) nhdV = 0 , where B{t)-.= {1 — t)Bi + tB 2 . (8.2) 

Then there exist a nonnegative integer n and a family {ujj(-)}f of continuous complex-valued functions 
on [0,1] such that for all t G [0,1], 


T,(B(t)) nV = {(.Vj(t)}f taking multiplicities into account. 


T 1 


see Section 3.1), the total multiplicity of quasi-eigenvalues 


Proof. Since F(--, B) ^ 0 for every B G 
of R in D is hnite. It follows from assumption (8.2), Proposition 4.2, the Weierstrass preparation 
theorem, and the Puiseux series theory (see e.g. |38l Theorem XII.2] and [HI [211 [23]) that for every 
s G [0,1] there exist £(s) > 0, a nonnegative integer and continuous in A := (s — £(s),s -f 
e:(s)) functions Q,sj{t), j = 1, ..., such that S(R(t)) fl D = for all f G A (taking 

multiplicities into account). 

It is easy to see that Ug does not depend on s. Let us denote this number n. Choosing a minimal 
hnite subcover from the open cover lJte[o 1] [O’ construct (using A and flsj(t)) 

a hnite partition 0 = to < ti < ■ ■ ■ < tk = I of the interval [0,1] and families {(-^i,j(-)}f=i, i = 0 , ..., 
k — 1, of continuous on [ti, A+i] functions with the property that T,(B(t)) HV = {ojij(t)}f^i for every 
t G [fi,fj+i] (taking multiplicities into account). Now one can continuously ’glue’ these families at ti, 
..., tfc-i to produce the desired family {(Uj(■)}". □ 


Proof of Theorem I g. 4 Let us take h > 0 such that the closed discs ©^(coo,^) are pairwise disjoint and 
all he in V. By Proposition |4.2[ choosing R> ||Ro||i; we can ensure that there exists a neighborhood 
W of Bq in the topology of the metric space (B/j,p) such that for all R G IT: 


(a) the total multiplicity of quasi-eigenvalues of R in each of the open discs ©^(cooj) coincides with 
the multiplicity of coqj for Rq, 

(b) the total multiplicity of quasi-eigenvalues of R in R coincides with that of Rq, 

(c) there no quasi-eigenvalues of R on each of the boundaries bdR, bdD5(a;oj), j = f) ■ ■ ■ 1 n. 
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B£(i?o) C W and take 61,2 satisfying (8.1). Then A C W. Denote 
are closed. By properties (a) and (b), S[A] flD = UA^. By (c), 


Let us take £ > 0 such that 
Aj := S[A] nD5(c(;oj)- Then Aj 
Aj C D5(ci;oj)- Proposition 8.3 and the fact that A is convex imply that each Aj is path-connected. 


This proves (l.i). 


Taking Uj G Aj such that Imca^ = max^^gA^ Imcj and Bj G A such that Uj G one can see 


that {uj, Bj) is a local minimizer for Dr. Since V is bounded, so are Aj. Hence, the existence of local 
maximizers for Dr can be shown in the same way. This proves (l.ii). The arguments for (l.iii) are 
similar. After all these considerations, statement (2) is obvious. □ 


9 Locally extremal uj are nonlinear eigenvalues 

Recall that the set of nonlinear eigenvalues was introduced in Section This section is devoted 
to the following theorem, which states that quasi-eigenvalues optimal in several senses lie in . 

Theorem 9.1. A number u is a nonlinear eigenvalue if at least one of the following conditions hold: 
(i) UJ is a guasi-eigenvalue of minimal decay for the freguency Recn, 

(a) Recj 7^ 0 and uj is a boundary point o/S[A], 

(Hi) there exists Bq such that {uj,Bq) is a local extremizer for Dr, 

(iv) there exists Bq such that (ca, Rq) is a local extremizer oflmuj for the freguency Reca. 

The proof given below combines the results of Sections [5](^ with the technique of local boundary 
points of S[A] introduced in [23] . 


9.1 Local boundary of E[A] and the proof of Theorem 9.1 


Definition 9.1 (cf. [23]). Let R G A and uj G S(R). Then: 

(i) the complex number uj is called a B-local boundary point of the image S[A] of A if a; is a boundary 
point of E[A n IT] for a certain open neighborhood IT of R in the norm topology of L^. 

(ii) UJ is called a B-local interior point S[A] if uj is an interior point of the set S[A fl IT] for any open 
neighborhood IT of R (in the norm topology of L^). 


This definition is a particular case of [23l Definition 4.1]. (In [23|, such uj are called strongly local 
boundary points.) Note that, if R G A, then a complex number uj G S(R) is a R-local interior point 
of S[A] exactly when it is not a R-local boundary point of S[A]. 

The following statement is a particular case of [23l Theorem 4.1 and Proposition 4.2]. 

Proposition 9.2 ([23|). Let Rq G A and uj G E(Ro). If cone [A — Rq] = C, then u is a 

Bo-local interior point o/S[A]. 


Consider the case u ^ iM. Combining Proposition 9.2 and Theorem 5.1 with the equivalence 


(5.15), one easily gets the following result. 


Theorem 9.3. //Rea; 7^ 0 and uj is a Bo-local boundary point o/E[A] for a certain Bo G A, then 
UJ G 


Now we are able to obtain Theorem 9.1| in the case Rea; = 0. Indeed, in this case, each of 
conditions (i)-(iv) implies that a; is a Ro-Iocal boundary point of S[A] for a certain Rq. Theorem 9.3 
complete the proof. 
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The proof of Theorem |9.1| for a; G iM follows from Proposition |7.1[ Indeed, it is enough to 
prove that (iv) implies u G (see Remark 6.1). Assume that (ca, Rq) is a local extremizer of Imca 


for the frequency 0. By Proposition 7.1, Bq coincides either with bi, or with 62 - 

Suppose Bq = 62 - Then y{x) = ^^9{x,u;b2) is a mode associated with the quasi-eigenvalue u 
of 62 - By (5.13) and Lemma(hi), 6 *^ > 0 everywhere on [ 01 , 02 ] except possibly one point. Hence 


y'^{x) G C+ a.e., and so 62 = 


bi + (62 - bi)Xc^ 


a.e. This means that y is an eigenfunction of 


(5.1), (1.2), (1.3) associated with the nonlinear eigenvalue u. 


If Bq = bi, we can put y{x) = e ™^^9{x,u;bi). This completes the proof. 


10 5-D reduction with the use of bang-bang equations 

10.1 Finding of optimizers from the nonlinear spectrum 

Assume that we have found the set of nonlinear eigenvalues. Then the quasi-eigenvalues of 
optimal decay can be calculated with the use of the following formula: 


/3min(«) = niin{—Imcj : a; G and Recj = a}. 


( 10 . 1 ) 


(the minimum exists for each a G ReS[A]). This statement follows from (5.3) and Theorem 9.1 As 
a by-product we get 

ReE[A] = ReE^'. (10.2) 


Roughly speaking, the nonlinear eigenproblem (5.1), (1.2), (1.3) excludes unknown optimal Bq 
(in any sense described above) from the optimization problem if corresponding optimal ujq is known. 
In this case, Bq can be recovered by one of eigenfunctions of the nonlinear problem. Rigorously this 
is formulated in the following result. 

Let us dehne a map y h->• B{y) from C[ai, 02 ] to Lj^(ai, 02 ) by 

B{y){-) = b,{-) + {b2{-)-b,{-))x^y{-)). 


Corollary 10.1. (a) Assume that statement (i) or statement (ii) of Theorem 9.1 hold. Then u is a 


quasi-eigenvalue of Bq ^ A if and only if Bq = B{y) a.e. on ( 01 , 02 ) for a certain eigenfunction y of 


(5.1), (1.2), (1.3) 


(b) If at least one of the statements (Hi), (iv) of Theorem 9.1 hold, then Bq = B{y) a.e. on ( 01 , 02 ) 
for a certain eigenfunction y of (5.1), 


(c) If Reu 7 ^ 0 and u is a BQ-local boundary point of S[A] for a certain Bq G A, then Bq = B{y) 
a.e. on ( 01 , 02 ) for a certain eigenfunction y of 0 , Q, Q. 


Proof, (c) By Proposition 


9.2 


cone 


dF{ujo',B()) 

dB 


[A — Bq] 7 ^ C. Combining this and Theorem 


5.1 


one 


obtains the desired statement. 

In the case Reco 7 ^ 0, (b) follows from (c). In the case Recn = 0, (b) easily follows from the 


arguments at the end of Section [9T| Note that in the latter case, Bq coincides either with bi, or with 
^ 2 - 

(a) The implication ’if’ is obvious. Let us prove the part ’only if’ in the case Re a; 7 ^ 0. Assume 
that CO is a quasi-eigenvalue of Bq G A. Since u G bdE[A], it is a Ro"local boundary point of S[A]. 
Thus, (c) implies the desired statement. In the case when Reco = 0 and co is of minimal decay, the 
implication ’only if’ follows from (b). □ 
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10.2 Nonlinear eigenvalues that are not of minimal decay 


Remark 10.1. If Bq is optimal in one of the senses of Corollary 10.1 then Bq is an extreme point of 


A. More precisely, Corollary 10.1 implies that after a possible correction of Bq on a set of measure 


zero, Bq{x) satishes at least one of the equalities Bo{x) = bi{x), Bq{x) = h 2 {x) at every x G ( 01 , 02 ). 
If, additionally, hi and 62 are piecewise constant with a finite number of intervals of constancy (see 
EZH for the precise formulation in 1-side open case), then Bq is so. The last statement follows from 
(For optical problems this means that Bq represents a multilayer resonator with a 


Lemma 5.2 


V 


finite number of layers.) 

Proposition 10.2. Let Rea; 7 ^ 0 and u G S(i?o) for a certain Bq & A that is not an extreme point 
of A. Then uj is an interior point of S[A]. 


The proposition follows from the combination of Proposition 9^ with Theorem 5T and its proof. 
Remark 10.2. Applying Proposition 10.2 in the settings of Theorems 8.1 and |8.2[ one can see that 


if UQ^j ^ iM, than the corresponding path connected component Aj (of the set S[A]) has a nonempty 
interior. 


8.2 


Example 10.3. Let Bq = b ^ b ^ S'lid letwoj be equal to 0Jj{b) of (3.7) for I < j < n. 
Let £ and 61^2 satisfy statement (1) of Theorem 
quasi-eigenvalue of minimal decay. 


Then there exist uj G 11“* \ iM that is not 


Indeed, by Remark 10.2, the interiors intAj of Aj are not empty. So for any a from their real 

and its proof imply the existence of a local minimizer {u},B) of 


projection Re[intAj], Theorem 8.2 


Imo; such that u) is not of minimal decay for a (see Fig. 1 in Section]^ 

10.3 Nonlinear eigenvalues as zeros of a function of three variables 


Theorem 10.4. Assume that bi ^2 satisfy (2.3) and 

hi{x) >0 a.e. on the set E (recall that E = {x : 61 (x) 7 ^ b 2 {x)}). 
Let a G [ 01 , 02 ] and cq,ci G C. Then there exists a unique solution y to the problem 

y{a) = Co, y'{a) = ci 

for the nonlinear equation 

y‘‘(x) = -z'-‘p{x) bt(x) + (fca) - 


(10.3) 


(10.4) 


(10.5) 


The proof includes a technical part that considers a number of cases for placing of co,i in C. It is 
postponed to Appendix |A.4[ 


Example 10.5 (Nonuniqueness). If assumption (10.3) is dropped, then the statement of Theorem 10.4 
is incorrect. To see this, consider the case when 612 are constants and 0 = 61 < 62 - Let G C_, 
o = Oi, Co > 0, and Ci = 0. Then, for small enough positive c, the functions y^ dehned by |/e(x) = Cq 


for X G [oi, 02 — s] and by |/e(x) = Co cos{zb 2 [x — 02 + e]) for x G [02 — s, 02 ] are solutions to (10.4), 


(10.5) on [oi, 02 ]. 


We assume (10.3) in the rest of this section. Assuming also ^ G ( 0 , 7 r], let us define the function 


0(x;^,z) as the solution to (10.5) satisfying the initial conditions 


0 (ai;^, 2 ;) = d,;0(ai;^, z) = -Rz/ic'^. 


( 10 . 6 ) 

























28 


Note that 0 satisfies (1.2) for any ^ (with u = z). Plugging 0 into the second boundary condition 


(1.3), let us introduce the complex-valued function 

19x0(02; 


Fniii^z) := 0 ( 02 ;^,^) + 


Zh'2 


for e C \ 0, 


3.1). 


and 0) := 1 - 1 - ^ (the logic of this definition is the same as in Section 

The following corollary essentially reduce the problem of finding of 
function /3min) to a finite-dimensional question of finding of zeros for a function of three real variables. 
Consider the fiber of the function F^\ over 0, 


(and so, of finding the 


F-i\0) := {(e, z) e (0, vr] X C : F„i(e, z) = 0}, 
and the projection of (0) to the z-plane 

pr^ F~^^{0) := {z e C : z) = 0 for certain ^ G (0, vr]}. 



without proof in 


11 A numerical experiment for constant side constraints 


To check applicability of the method of Section [10.3 to calculation of quasi-eigenvalues of minimal 
decay, we take ui = 1, 1^2 = +00, Oi = — 1 , 02 = 0 , and the constant functions bi = 90, 62 = HO. 
Then Fni(^, z) = 0(0; z) and = {z E C_ : Fni(^, z) = 0 for certain ^ E (0, tt]}. 

To find numerically, we take small £ and approximate the sub-level set ^ni(£) '■= {z E C : 
inf^g(o, 7 r] |-^ni(OI — Oy its discrete versions Z{hR,hi^ha.rg,e) consisting of numbers z = zi F ' 1 Z 2 
with Zi E hjiZ, Z 2 E hjZ, such that maxi<„<Ar |Tlii(nharg, ^)| < £. The angular step harg is connected 
with N by harg = tt/N. 

The value of F^i{^, z) = 0(0; z) is computed by the shooting method applied to the solution 0 
of the equation 

y"(x) = -z^y{x) bi + [62 - bi]x^^iy^ix)) . ( 11 . 1 ) 


This does not require finite-difference approximation. Indeed, for G C_, the equation (11.1) has 
constant coefficients on each interval where 0 ^ stays in one of half-planes C± and admit an explicit 
analytic solution depending only on the initial values in the left end of this interval. These values 
and the lengths of intervals of constancy can be computed iteratively starting from initial conditions 


( 10 . 6 ). 


The result of such a computation of Zni{e) with the level e = 5x 10“® in the domain V = {Re z G 
( 0 , 1 . 2 ), Im^ G (— 0 . 015 , 0 )} are plotted in Figure|^(a). We take N = 360 . The other step parameters 
are changed through V from {hji, hj) = ( 2.4 x 10 ^, 5.5 x 10 “®) to {Hr, hi) = ( 2.7 x 10 “®, 2.3 x 10 “"^) 
to provide better resolution for the regions where Z^\{e) is concentrated and to reduce computation 
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time in the regions where the values of function Ai{z) := maxi<„< 7 v |-Fni(’^harg, ^)| are essentially 
greater than e. 

Comparing the shape of this approximation of fl V with the fact that bdS[A] \ iM C 
(see Theorem 9.1 (ii)) and Theorems 8.18.2, we conclude that the considered example hts into 
the description of the small contrast case at least in the domain V, and that Z.^\{e) provides an 
approximate shape not only for fl but also for bd S[A] fl V. The form of bd S[A] fl T> together 
with Theorem 8.2 makes it possible to find approximately the part E[A] flP of the set of admissible 
resonances, which consists of four quasi-eigenvalue ’clouds’ Ai, ..., A 4 (see Remark 6.2). 



Figure 2: MATLAB computation of: (a) E"' in the rectangle V, (b) the graph of M.{z) := maxi<„<jv |Aii(n/iarg, z)\ 
in the neighborhood I ?4 = {Rez G (1,1.2), Imz G (—0.025,0)} of A 4 . 


While the shape of in the rectangular neighborhood = {Rez G (1,1.2), Imz G 

(—0.025,0)} still resembles a boundary of a closed bounded path-connected component A 4 of S[A] 
(and even suggests that A 4 is simple connected), the graph of function M.{z) exhibits a close to peri¬ 
odic pattern of ridges and canyon in the regions with Im^; close to —0.007 and —0.013, see Figure 
(b). For the study of A4, the value of harg = vr/360 is too large. When harg is decreased to ^ x 10“^ 
while keeping > 10“^ and hj > 10“®, only one sharp canyon that corresponds to bdS[A] fl 
remains. 

The values of the minimal decay rate function /?min for « in the intervals Re Ai, ..., Re A 4 can 
be computed using (10.1). We have done these computations for some of the values of a, see Table 
1. The step hj is 5 x 10“^ for the first and second values of a and hj = 5 x 10“® for the others, h, 


varies from x 10 ^ to 3 !^ x 10 


arg 


a 

0.14977 

0.14983 

0.16557 

0.44931 

0.46050 

0.49674 

/3min (*^) 

0.009119 

0.009120 

0.01115 

0.00910 

0.00846 

0.01115 


a 

0.74884 

0.77200 

0.82790 

1.04838 

1.08800 

1.15905 

/^min (*^) 

0.00910 

0.00766 

0.01113 

0.00909 

0.00689 

0.01109 


Table 1: Some of values of the minimal decay rate function 0 m\-n (<y)- 


For a = 1.088, we have found the structure B generating the quasi-eigenvalue of minimal decay 
a — i/d(Q;). It consists of 7 layers. With xq = —1 and X 7 = 0 for the endpoints of the interval 
[oi, 02 ] = [—1, 0], the computed intervals of constancy of B can be written in terms of the boundaries 
of the layers (switch points) Xj, j = 1,..., 6. Namely, B{x) = 110 for x G {x 2 k)X 2 k+i)) A: = 0,..., 3, 
and B{x) = 90 for x G {x 2 k+i^ X 2 k+ 2 ) for k = 0,..., 2. Approximate values Xj for Xj are given in 
Table 2 together with absolute error estimates Sj. 
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J 

1 

2 

3 

4 

5 

6 

Xj 

-0.86237855 

-0.71669445 

-0.57371202 

-0.43130025 

-0.28563822 

-0.15697427 


0.00000028 

0.00004800 

0.00000900 

0.00007000 

0.00002000 

0.00011200 


Table 2; Switch points and their absolute error estimates for B{x) of minimal decay for a = 1.088. 


One of possible improvements of the proposed method may be concerned with calculation of zeros 
of the function Indeed, we used the simplest straightforward approach based on sub-level 

sets of the function Ai. More elaborated approaches that do not neglects the information contained 
in a.Tg F^\{nharg, z) should be more efficient. 


12 Discussion 

Earlier Engineering and Numerical Optimization papers were concerned mainly with the case when 
the constraint functions bi and 62 are constants. In this section, we (sometimes speculatively) compare 
our results for this case with earlier suggestions for optimal designs and with discussions of the proper 
statement of the problem. 


12.1 Alternating almost periodic structures with a centered defect. 

Contemporary designs of high-Q optical cavities usually involve incorporation of defects into a peri¬ 
odic structure composed of two materials [h [33113S] • For cavities with 1-D geometry [1111331 ESI E], 
the base periodic structure consists of alternating layers (say with widths ii and ^2 ) of two materials 
with different permittivities (say hi and 62, respectively). The Physics justihcation for such designs 
is that the periodic structure forms a distributed Bragg reflector with high reflectivity for certain 
bands of frequencies (stopbands). When a defect is introduced, the waves with such frequencies are 
expected to accumulate and to be well confined in the region surrounding the defect m Chapter 
4], [32]. If the fabrication process involves material with permittivities from hi to 62, Remark 10.1 


implies that the best confinement has to be produced by structures B that consist only of layers 
with B{x) = hi and B{x) = 62 (for one-side open case, see [211122] )• In other words, both the base 
structure and the defect have to consist of the two materials with extreme allowed permittivities. 

Basing on various experimental and computational approaches, two types of designs with defects 
were proposed: 


(a) a localized defect in a center of a periodic structure, e.g., [DEUi, 

(b) the size-modulated 1-D stack [331 ESI 13] • 

Designs of type (b) assumes gradually changing deviations from periodicity in the widths of layers. 

For the structure of minimal decay computed in Section]^ (see Table 2), one can observe that 
it is close to alternating periodic with a period consisting of one layer with B{x) = hi and the other 
layer with B{x) = 62 (since the structure is finite, either the first, or the last layer has to be treated 
separately). The widths of layers [xi^xfl, [0:3, xj, [x^,xf\ with B = hi gradually decrease from 0.146 
to 0.129, and the widths of layers [x 2 k, a^2fc+i], fc = 0,..., 3, corresponding to R = 62 gradually increase 
from 0.138 to 0.157. This resembles designs of type (b). 

Let us note that the 1-D structures considered in [l9l[331EniE] are symmetric. If we consider 
the numerical experiment of Section in the context of optimization of odd modes of a symmetric 
structure in the interval (—1,1) (cf. Example |6.4[ ), then the structure B has to be extended to (0,1) 
by B{±x) = B{x), X G (—1, 0). Then the layer [xe, —Xg] is approximately twice larger of the nearby 
layers and so the structure does not resemble the size-modulated 1-D stack designs of [331 EHl E] . 
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The symmetric structure B can be considered as a combination of designs (a) and (b). It has a 
localized defect consisting of the central layer [xq,—x^, and this defect is surrounded by alternating 
layers with gradually changing widths. 

In [321 [2ni E] , the quadratic deviation laws were employed to vary the widths of layers. From the 
point of view of analytic optimization, it would be interesting to hnd dependence of the base periodic 
structure, the parameters of the defect, and width deviations on the frequency of corresponding 
optimal resonance. 


12.2 Global minimization of decay rate. 


Well-posedness of global maximization for Q-factor related functionals (without frequency restric¬ 
tions) was discussed in [191129]. For 1-D resonators with constant side constraints and hxed length 
02 — oi, the global maximizer for Q-factor does not e xist . Indeed, if we take a homogeneous structure 
B{x) = b (with constant b G [&i,& 2 ], b ^ ^^ 12 ); then (3.7) implies that Q{uJn) = go to -|-cxd as 

|n| —)■ 00 . 

Below, basing our intuition on the numerical experiment of Section H, we discuss the question 
of global minimization of the decay rate | Ima;| (cf. [T5]i. 

It can be seen from Figure]^ (a) that quasi-eigenvalue ‘clouds' corresponding to higher frequencies 
are ‘wider’. It is natural to expect that for | Ima;| great enough they intersect each other and form 
one large unbounded cloud. This conjecture is supported also by our study of clouds’ projections 
(see Section 6.2 and Appendix |A.3[ ) and by the results on existence of multiple quasi-eigenvalues (see 
Section 3.2). In particular, we already know that the minimal decay function is dehned for large 
enough |a| (Proposition |6.5| (hi)). 

It also can be seen from Figure (a) and Table 1 that each of subsequent clouds Ai, . .., A4, 
comes closer to the real line than previous. (The values /3(0.14977) ~ 0.009119, /9(0.4605) ~ 0.00846, 
/3(0.772) 0.00766, and /5(1.088) ~ 0.00689 approximately corresponds to minima of decay rate for 

the clouds Ai, ..., A4.) Basing on the above observations, we propose for the case when 61,2 are 
constants and 0 < 61 < 62 < +00 the following 


Conjecture 1. The inequality liminf|Q|_,.oo/3 (q!) < /9(ao) holds for all admissible frequencies uq. (In 
particular, the global inhmum of the decay rate inf^gs[A] I Inia;| is not achieved over A.) 

A stronger form of this conjecture that assumes existence of hm|Q,|^oo (d{a) could help to explain 
the sliding effects in gradient ascent numerical experiments of [151 [T9|. Namely, in these experiments, 
iteratively improved resonances a;!"'] were approaching to the real axis with simultaneous growth of 
frequency Reca^^l This could happen if for some iteration uq reach (or come very close to) 
the Pareto optimal frontier {a — i/3(a) : a G ReE[A]} and then move along this frontier to 00 
with growing Reca and decreasing | Ima;|. In our opinion, the comparison of results of simulations 
[31 Ha UHl E3] suggests that lim|Q,|^oo/d(«) = 0 . 

It is also interesting to consider the value of liminf|Q,|^oo/d(Q^) from the point of view of the 
question of uniform separation of the set of resonances E(R) from the real line (the question of 
exponential energy decay in the sense of [SllTlEI). The absence of quasi-eigenvalues in the certain 
strip {—/do < Imz < 0 } is presently known only under additional conditions involving bounds on 
total variation of R [7]. So it is difficult to expect that, for certain (3q > 0, the set S[A] does not 
intersects with the strip {—/do < Im^ < 0} in the case of the admissible family A restricted only by 
constant side constraints. 
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A Appendix 


A.l Analyticity, the proof of Lemma 3.2 


(i) follows from (ii) and ( 3.4| ). Let us prove (ii) for {z,B) i—Consider the Maclanrin 
series (in e C and B G Lj.) 


'h(a;, 2 ;; B) = {x - oi) - B)z^ + B)z^ - ip'i{x] B)z^ + ..., 

pX 

'4)q{x-, B) = X - ai, ipj{x]B)= {x - s)ipj_i{s]B) B{s)ds, jeN, 


fai 


as a C[ai, a2]-valned series. One can see that 

\^jJj{x;B)\ < il)j{x] \B\) < il)j{a2] |i?|) for all x G [01,02]. 

It can be shown in the way similar to [HI Sect. 2 ] that 

w 2^(02 -Oiy+i5||{ , 

i’j{a 2 ] \B\) < --, j G N. 


(A.l) 

(A.2) 

(A.3) 

(A.4) 


Hence, the series (A.l) converge nniformly on every bonnded set of C x L^. So (A.l) dehnes an 


analytic map from C x to (^[O,^]. Since (A.l) satisfies dly = —z'^By, y{ai) = 0, and dxy{ai) = 1, 
we see that d/{x,z;B) = 'ip{x, z] B). The eqnality 5^-0 = —z'^B'ip yields that G and, moreover, 
that {z,B) I—)■ z] B) is an analytic map from C x to L^. Hence, (A.l) defines an analytic 

map from C x to This map is also bonnded-to-bonnded dne to (A.3), (A.4) and the 

eqnalities 9^-0 = —z'^B'ip, 'ip{ai) = 1 . 

Similarly one can prove that {z, B) 1 —)■ 0(-, z] B) is an analytic and bonnded-to-bonnded map from 
C X to (cf. [iHl |23]). Finally, since 6 = ip — izui'ijj, we obtain the same statement for the 
map {z,B) I—)■ 9{-,z;B). This completes the proof. 


A.2 Directional derivatives, the proof of Lemma 4.3 


Differentiation w.r.t. B. For any z G C, the fnnctions 9(x, z] B) and dx9{x, z] B) satisfy (3.2) and 


dx9{x,z;B) = —izvi — z^ / B{s) 9{s^z]B) ds. 


(A.5) 


' 0.1 


To hnd directional derivatives §§(V) and we differentiate these eqnalities by dehnition 

using Lemma [3^ We get 


d9{x, z] B) 


dB 


{V) = —z^ {x — s)B{s) 


'ai 


d9{s, z; B) 

an 


{V)ds — z'^ / {x — s)V{s)9{s,z;B)ds, (A.6) 


'01 




dB 


dB 


It follows from (A.6) that iY) is the solution y{x) of the initial value problem 

y"{x) + z‘^B{x)y{x) = /(x), y{ai) = 0, y\ai) = 0 


(A.7) 


(A.8) 


with /(x) = —z‘^V{x)9{x,z;B). Further, (A.7) can be rewritten as = y\x) (as a 

by-product, we get ^ d^e(x,z,B) ^ ae(x,z,B) Solving (A.8) by variation of parameters, one 
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can find y,y', and, in tnrn, y{x) + For 2 ; 7 ^ 0, [y{x) + eqnals 


0(02) H- dx 0 {a 2 ) 

ZV2 


f{s) ij{s,z]B) ds + 


'ai 


^(02) H- dx'ip{a 2 ) 

Zl>2 


Snbstitnting /, and z = uj E S(i?) , we get 


f{s) 9{s,z]B) ds. 


fai 


9{a2,uj;B) + - dx9{a2,uj; B) = 0 

UJl'2 


and — (V) = -uj 


^(02) H- d^'ilj{a 2 ) 

ZV 2 


To obtain (4.3), it remains to note that 


'■02 


9‘^{s,u:-B) V{s) ds. 


'ai 


^( 02 ) +- d^'4}{a2) = - ^ 

ZV2 Oa;9{a2) 


The last identity easily follows from the constancy of the Wronskian 


9 ip 
dx9 dxp) 


= 1 . 


Differentiating (3.2) and (A.5) w.r.t. z, we see that 8^9 is given by the solntion y of the 
problem 

y''{x) + z^B{x)y{x) = f{x), y{ai) = 0, y\ai) = -ivi 
with f{x) = —2zB{x)9{x, z; B) and that dzdx9 = 8x8^9. Hence, for z = a; G S(H), 


8zF{uj;B) = 8z9{a2,uj;B) + 


i8z8x9(a2,uj;B) i8x9{a2,u-, B) 


00 U2 


Uj‘^l'2 


= y{a2) + 


iy'{a2) i8x9{a2,u:]B) 


UIZ 2 


U‘^U2 


Solving for y by variation of parameters as before, one can easily get (4.2). 


A.3 Admissible high frequencies, the proof of Proposition 6.5 (iii) 


It is enongh to prove the statement of Proposition 6.5 (iii) nnder the assnmption that 61,2 are constants 
satisfying hi < 62 - We consider in details several essentially different cases of arrangement of the 
intervals [z/i, z/ 2 ] and [5}^^, and very briefly the other cases (for which the proof can be obtained 
by simple modihcation of argnments). 

Step 1. For the cases when 0 < 9^^ < < ui or z /2 < 9^^, let us show that 


ReS[A] D < |a| > 


TT 


(02 — fll)(^2 


1/2 


K 


1/2 


- bT 


(A.9) 


Indeed, for each n G N consider the intervals Sn covered by Rea;„( 6 ) when a constant b changes in 
the range (&i, 62 ]- By (3.7), Sn = a^-ai h the right end of Sn is > than the left end 


of for all n > no G N, then ReS[A] contains the set = 


obtains (A.9). 


nriQ 


(«2- 01)^2 


YJ 21 + 00 ^ and also its 


reflection w.r.t. 0, see Lemma 3.3 (ii). Estimating no from the inequality hi ^'^^n > 63 + 1), 


one 
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Step 2. Slightly modifying arguments of Step 1, one can prove that: 


ReS[A] D < |a| > 


ReS[A] D < |a| > 


TT 


(02 - 


1 

+ 




1/2 


- hT 2 


TT 


(02 - 01)62'^^ ' 2 


1 

+ 


K 


1/2 


1 

+ 


- hT 2 


ReS[A] D < |a| > 


TT 


(02 - 01)62'^^ 

and that, 

in the cases 0 < 61 < 63 ^^ = and vi < h\^‘^ < = z /2 


if vi < b\^'^ < buj"^ < 1 ^ 2 ; (A. 10 ) 

if Pi = < b^"^ = 1 / 2 ; (A. 11 ) 

if 0 = b\^‘^ < b^^ < Pi] (A. 12 ) 


(A.9) and (A. 10 ), resp., are valid after the change of |a| > ... to |a| > .... (A.13) 

1/2 1/2 

Step 3. Let us show that in the case Pi < b^ < P 2 < b^ , 


Re E[A] D < |q:| > 


TT 


\br+br 


b'p - ip 


(02 - 01 ) 62 '^^ 

Now, the sets Sn covered by Rea;„(fe) have gaps. Indeed, assume that bi > 0. Then 

71 [ n + 1/2 n + 1/2 


(A.14) 


TT 


Sn — , , 

(02 - Oi) 


n n . 
^,7- U 


bV^'^2 (02-ai) 


^2 bl 
-1/2 


1/2 


The gap is closed by the sets iS„_i an d cSn+ i as soon as (n — l/2)6A/^> {n + 1)62 The latter 


1 ^ 1/2 ^ 1/2 
2^2 ^^1 


_^ ~p(-/ 1 /2 

is valid if n > This proves (A.14) for bi 7 ^ 0. To prove (A.14) in the case D = b-/ , it is 

+ — bl 

enough to take the limit 61 —)■ 0 . 


Step 4- Modifying arguments of Step 3, it can be shown that 

1 


ReS[A] D |a| > 


TT 


(02 - 01)62'^^ ' 2 


+ 


351/2 

2^1 


Re S[A] D < |a| > 


TT 


(02 - ai)bl^‘^ 


bl'^ - bT 


ReE|A]3pa|> 

((2-2 ^ 1)^2 


ReS[A] D s |a| > 


27r 


(02 - ai)&2'^^ 


if 0 < 6 ^'^^ <Pi< by^ < P 2 ] (A. 15) 

if 0 < by < Pi < P 2 < (A. 16) 

if 0 = <Pi < by < P 2 ] (A. 17) 

if 0 = by < Pi < P 2 < by. (A. 18) 


Thus, (6.2) follows from the above inclusions. 


A.4 Solvability of bang-bang equations, the proof of Theorem |10.4 

It is enough to consider the case z'^ G C_ (see the explanations to (5.2) in Section]^. 


Lemma A.l. Let (10.3) be fulfilled. Let G C_. Then there exists e = e{a) > 0 such that problem 


(10.5), (10. 4 ) has a unique solution y on {a — e, a + e) H [oi, 02 ]. 
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Proof. Firstly, if hi = 62 in a neighborhood of a, the statement is obvious. 

Next, to be specihc, assume that a < 02 and that for any 5 > 0 the measure of the set of 
X G {a,a + 6 ) such that 61 (x) 7 ^ 62 ( 2 ^) is positive. Then (10.3) yields 

for any h > 0, measjx G (a, a + h) : bj{x) > 0} > 0, j = 1, 2. (A. 19) 

It is enough to prove the lemma for the interval [a, a + e) with small enough e > 0. 

Uniqueness. Assume that there exists at least one solution y to (10.4), (10.5). Our goal is to 
prove that there exist e > 0 and Bq{-) such that for any solution y to (10.4), (10.5) the function 
B{y){-) coincides (in L^-sense) with i?o(') on (a, a + e). This and the structure of equation (10.5) 
immediately implies uniqueness. 

Case 1. If Co = Cl = 0, any solution y to (5.1), (10.4) is trivial and Biy) = hi. 

Case 2. If Cg ^ M or d^y'^^a) = 2coCi ^ M, then, for a certain nonempty interval (a, a + c), the 

sign of Im|/^(x) is completely determined by cuandci. This easily gives the desired statement with 
Bq equal to one of the constraints (either hi, or 62 ) on (a, a + e). 

Case 3. The case when cq = 0 and cf ^ M is analogous to the Case 2 . 

Case 4 . Assume that Cq G M \ {0}, CqCi G M, and G M. Then either cq, Ci G M, or cq, Ci G iM. 
Moreover, any solution y stays on the same line (M or iM) since G M. Therefore y‘^{x) G M and 
B{y) = hi. 

Case 5. The case when Cg = 0, c^ G M \ {0}, and G M is analogous to Case 4. 

Case 6 . Assume that Cg G M \ {0}, cqCi G M, and z"^ G C_. Consider an arbitrary continuous in a 

vicinity of a branch arg^ y{x) of oxgy. Since cqCi G M, the arguments of the proof of Lemma 5.2 (i), 
(iv) imply 


dx arg^ y = 


\m{yd:,y) 

\y? 


and \m.[y{x)dxy{x)] = —Imz^ / |?/(s)ff3(|/)(s)ds. 


It follows from these equalities and (A.19) that arg^ j/(x) > 0 on (a, a + Si{y)) with a certain 
Si{y) > 0. So for any solution y, 


ini;: 


on (a, a + Ci(?/)). 


(A.20) 


This easily implies that it is possible to choose ei{y) in (A.20) independently of a solution y. 
Case 7. Assume that cq = 0, ci G M \ {0}, and z"^ G C_. Since 


y'{x) = Cl - / ?/(s)i3(?/)(s)ds = ci - z"^ 



B{y){s)y\t)dtds, 


(A.21) 


we see that 


lmy\x) = 



B{y){s) lm[z^y\t)]dtds. 


(A.22) 


Hence, for x slightly greater than a, the complex number y'{x) lies in the same open quadrant as 
Cl —i Im[z^ci]. Due to the equality y{x) = y'(s)ds, the numbers y(x) he in the same open quadrant. 
So the sign of Ci completely determines B{y){x) on (a, a + ei{y)) with a certain ei{y) > 0. As before, 
it is easy to see that Si{y) can be chosen independently of a solution y. 

Case 8. When cq = 0, ci G iM \ {0}, and z"^ G C_, one can use arguments similar to Case 7 with 
the change of (A.22) to 


Rej/'(x) = - 

This completes the proof of uniqueness. 



B{y){s) Re[ 2 :^|/'(f)]dMs. 
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Existence. The above arguments also provide an algorithm for assigning Bq{x) the values of 
bi{x), or b 2 {x) for X > a such that there exists a solution y to y” = —z’^B^y on a certain interval 


(a, a + e) satisfying the initial conditions (10.4) and additionally Bq{x) = B{y){x) on (a, a + e). Then 
y is a. solution to (5.1), (10.4) on (a, a + e). This completes the proof. □ 


Assume that (5.1), (10.4) has a unique solution y on {xi,X 2 ) with xi < a < X 2 < +oo. 

Lemma A.2. There exist ^i'^x^x 2 -u'i.^) ■ 

Proof. The statement follows from the Bellman-Gronwall lemma, from the integral reformulation of 


problem (10.5), (10.4) given by the equation 


y{x) = Co + ci(x - a) - / (x - s) y{s)B{y){s)ds, 


(A.23) 


and from the fact that \B{y){-)\ < b 2 {-) G L4(a;i,X2) for arbitrary y. 


□ 


This lemma, the equality y\x) = ci — y{s)B{ii){s)ds, and (A.23) imply that y and y' can 
be considered as absolutely continuous functions on [ 0 , 0 : 2 ]. This arguments allow us to extend y on 


[0,02]. The resulting solution is unique due to Lemma A.l This completes the proof of Theorem 

MM 
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